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ABSTRACT 


This thesis is concerned with elastic critical loads of single story, 
Single bay, rectangular rigid frames subjected to uniform vertical load and 
concentrated horizontal load. These frames are hence subjected to axial loads 
and primary bending moments. The frames investigated have hinged bases, fixed 
bases, or base fixity intermediate between these two extremes, Slope-deflect- 
ion equations, modified for the effects of axial loads, are used for the stru- 
ctural analyses of the frames and an iterative technique involving the use of 
an electronic digital computer is employed to solve the governing equations 


for each frame, 


Critical loads are defined as the highest load on a load-deflection 
plot for the frame. The results indicate that the fixed-base frame will sus- 
tain an elastic critical load about four times greater than the correspond- 
ing hinged-base frame, for a constant value of horizontal load, When horiz- 
ontal load is present, both frames exhibit large sway deflections before the 
critical load is reached. For both frames, the elastic critical load decreases 
with increasing horizontal load, the frame dimensions and ratio of member 
rigidities remaining constant. The elastic critical load also decreases as 
the length of beam increases with respect to the length of columns, the 
horizontal loa’ and ratio of member rigidities remaining constant. Further- 
more, the elastic critical load increases as the beam rigidity increases 
with respect to the column rigidity, the horizontal load and frame dimensions 
remaining constant. The presence of even a small amount of base restraint 


increases the critical load appreciably above that of a hinged base frame, 
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This amount of base restraint is available under presently used column bases 


that are commonly assumed in design to be hinged bases, 


This method of solution may be readily adapted to stability analysis 
of other practical building frames subjected to axial loads and primary bend- 


ing moments, 
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NOMENCLATURE 


C = Non-dimensional carry-over factor 
E = Modulus of elasticity of frame material 
EI = flexural rigidity of frame member 


Flexural stiffness of coluwm members 

Flexural stiffness of beam 

Flexural stiffnessof restraining beam 

FEM = Fixed end moment 

I = Moment of inertia. of cross-section of frame member 
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: = ate factor of frame member = SEI/L 


= He of frame member 

L, - Length of column members 

L, - Length of beam 

L,,~ Length of restraining beam 

My = = Moment at end + of member lj 

P = Axial load 

= Simple Euler load of frame member 
(P, ), = Simple Euler load of colum members 
(P.), = Simple Euler load of beam 

RP = ®/,= Ratio of horizontal to vertical load 


REI = El,/EI, 


S = Non-dimensional stiffness coefficient 
w = Uniform lateral load 

x = Coordinate direction 

y = Coordinate direction 

0 = pee rotation of frame member 

0; = Rotation of joint Ut 


<= Total vertical load/(P, ), =wL,/(P, »; 
® = Total horizontal load/(P,), = RP x 


A= L/L, 
¥ = (EI/L), /(EI/L), 
N= (EI/L),, /(EI/L), = El,,/El, 
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CHAPTER I 
INTRODUCTION 


1.1 STRUCTURAL DESIGN CRITERIA’ 

The use of limit design in present-day engineering practice req- 
uires the satisfying of three main criteria for structural damage. One of 
these is the ultimate or plastic strength criterion, which is satisfied by 
insuring that the maximum load on a structure times a safety factor does 
not exceed the ultimate or plastic strength of the structure. A second criter- 
ion is that the buckling strength must not be exceeded, which has led to the 
concept of design for stability. Buckling occurs in the elastic or inelastic 
ranges when the structure passes from one deformation configuration into an- 
other without a change in the load. A third criteria to be satisfied is that 
the deformation of the structure must not exceed specified limits, in order 


that the serviceability of the structure is not impaired. 


The usual design procedure is to proportion the members on the basis 
of their ultimate strength, and then perform separate checks to insure satis- 
faction of the stability and deflection criteria. The individual members 


and the structure as a unit must satisfy all criteria. 


1.2 CONCEPT OF FRAME INSTABILITY 
Structural design for stability requires consideration of three general 
types of instability, these being: 


(i) Local buckling of elements of the cross-section of individual 
members, 
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(ii) Instability of individual members under the action of axial force 
and bending moment. This category includes instability of beam - 
columns due to excessive bending, and lateral-torsional buckling. 

(iii) Overall instability of the entire structure, often referred to 

as frame instability. This arises from unfavorable deflections 
in the frame as a unit, the mutual interaction of all members 
being the important feature. 

The first two types of instability have been extensively studied for 
many years, and satisfactory design recommendations have evolved from this 
work. The subject of frame instability has become more prominent with the 
advent of plastic design. Limitations on the application of the plastic the- 
ory have arisen from the difficulty of dealing with instability problems in 
partially plastic structures. Furthermore, modern architectural requirements 
have resulted in the use of frames virtually unbraced by wails or partitions, 
with light floors, and non-structural fireproofing. Without the built-in 


safety factors of conventional construction, frame instability becomes an 


important consideration. 


Although frame instability is sometimes of importance even in the 
elastic range, it is frequently ignored in design. The complexity of the 
problem may be the main reason for the comparatively small amount of atten- 
tion devoted to it, and this complexity increases greatly when plastic deform- 


ation occurs. 


Frame instability can be illustrated as in Figure 1-1 for a portal 
frame which is not prevented from sidesway. Two types of failure may occur: 

(A) When a symmetrical frame is loaded by symmetrical vertical forces, 
it is possible that the frame may pass from a symmetrical stable deformation 
configuration to an unsymmetrical, unstable configuration (cases 1 and 2). 


At this instant, the total resistance to any lateral force or lateral move- 
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ment becomes zero, The load-deflection curve is characterized by a sudden 
shift from a situation where the applied load can increase with little or no 
deflection to one where large deflections develop without an increase of 
load. The behavior is analogous to that of a centrally loaded column, in 
which bifurcation of the equilibrium position is possible at a certain criti- 
cal load. At a point of bifurcation, it is theoretically possible for either 
the stable or unstable configuration of the frame to be in equilibrium with 


the critical load. 


It should be noted that if the frame is prevented from sidesway, 
the critical load would occur when the frame transforms to a symmetrical 
unstable configuration. This load is always higher than the sidesway critical 


load, and hence does not govern for frames free to sway. 


There is a distinct difference in the nature of cases 1 and 2. In 
case 1, where the loads are applied directly to the columns, the frame the- 
oretically carries no primary bending moment before buckling. Therefore, 


only the column action alone needs to be considered in the buckling analysis, 
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However, most rigid frames are primarily designed to support loads such as 
in case 2, where all the members are subjected to both axial force and bend- 


ing moment before buckling. 


(B) When a frame (symmetrical or unsymmetrical) is subjected to a 
combination of horizontal and vertical forces (case 3), or when an unsym- 
metrical frame carries beam loads, the frame deforms laterally upon the app- 
lication of the first load. The resulting change in geometry may alter the 
carrying capacity of the individual columns, since the colum top is no longer 
diréctly over the column base and hence additional bending moments are intro- 
duced by the vertical load. The whole structure becomes unstable in this 
deformed position much like an eccentrically loaded colum. At the critical 
load, the structure continues to deform with no increase in load. Very little 


information is available concerning this type of failure. 


The above discussion has inherently pointed out the three main causes 


of sway collapse. These are: 
(i) vertical loads taken directly on columns. 


(ii) vertical loads on beams causihg primary bending. 
(iii) horizontal loads (wind loads). 


The first item (i) is the most important for tall building frames, 
since large axial forces cause a drastic reduction in the stiffness of the 
columns, The last two items are usually less important, although the exact 


extent of their influence has not been established, 


The case of a frame subjected to a combination of horizontal and verti- 
cal loads as in case 3 represents a situation where all three causes of sway 
collapse are present, The critical load occurs when the horizontal deflect- 


ion of the column tops increases with no increase in load. A linear elastic 
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analysis of the structural response of this frame does not exhibit such a 
maximum since actual collapse is due to factors not included in this type 
of analysis. These include primarily the effect of axial load on the stiff- 
ness of the frame members, the effect of plastic behavior, and the effect of 


finite displacements of the frame. 


1.3 RELATION BETWEEN ELASTIC AND PLASTIC ANALYSES 

In the calcuation of the elastic critical load, i.e, the critical 
load considering the frame to remain completely elastic, the effect of axial 
load on stiffness and the effect of finite displacements can be considered in 
the analysis. The calculation of the simple plastic collapse load by the 
usual simple plastic theory considers only the property of plasticity, although 
the effect of finite displacements can be considered by an extension of the 
simple plastic theory. (1) The relation of elastic and rigid plastic behavior 
to elastoplastic behavior can be most adequately described by reference to 


Figure 1-2, 


Load P 
Elastic Critical 
Load (P,, ) 


Perfectly Elastic 
=> Collapse 


Simple Plastic Collapse 


Rigid Plastic Collapse 
Load Including Effect o: 


Elasto-Plastic Critical Load pete t+ called 0 
Finite Displacements 


= True Failure Load (P, ) 


Elasto-Plastic 
Collapse 


Deflection Parameter 


FIGURE l<2 
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A perfectly elastic analysis considering the effect of axial load 

on stiffness and the effect of finite deformations produces a load deflection 
curve as shown. The elastic critical load P, then is the highest point on 
this load-deflection curve. The simple plastic theory yields the plastic 
collapse load (P, ) which is constant and independent of deflection. If the 
effect of finite displacements is included in the plastic analysis, a plastic 
collapse load is obtained which depends upon the deflection parameter. The 
true failure load is thus lower than that predicted by the simple plastic 
theory, and occurs due to instability of the elasto-plastic structure before 
a complete mechanism has been formed, The exact shape and position of the 
true (elasto-plastic) collapse curve is difficult to assess, but has been 


discussed by Horne. (1) 


To obtain an insight into the reason for an instability collapse 
before a mechanism has been formed, Wood(2) has proposed the concept of 
a "deteriorated elastic structure", From an elastic stability point of 
view, the plastic zones in the structure are incapable of developing add- 
itional resistance to displacement. The deteriorated elastic structure 
at collapse can therefore be obtained by removing all yielded material from 
the frame and replacing plastic hinges by frictionless hinges. The elastic 
critical load of this deteriorated structure is called the "deteriorated 
critical load". A frame analyzed by Wood showed that the deteriorated critical 
load drops drastically as successive hinges are developed. Presumably, a 
reasonable estimate of the true failure load of the structure would be the 


“Ndeteriorated critical load" of the structure for which all but the last 


hinge has formed. 
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The concept of the deteriorated critical load clarifies thinking 
on the subject of elastic-plastic instability of structures, However, 
it is of limited assistance only since the actual calculations are very 
difficult and time-consuming. For this reason, empirical approaches have 
been attempted to relate the true failure load (Pp ) to the elastic critical 
load (P- ) and the simple plastic collapse load Ger Merchant (3) has sug- 


gested a Rankine type formula as in Figure 1-3, 


FIGURE 1-3 


& 


" Rr 


Experimental evidence to date has indicated that this formula is 
most successful when P,/P,, is small and the true collapse load is close 
to P,. For larger values of P, /P,, the formula is conservative, It is 
expected then that this formula is most applicable to low buildings, where 
in fact the elastic critical load is usually much higher than the collapse 
load predicted by the simple plastic theory. For taller buildings, the for- 
mula would provide a conservative estimate or lower limit to the true fail- 


ure load, 


1.4 LITERATURE REVIEW 

Many investigations into the stability of frames have been carried 
out, almost all of this work being done under the assumption that the frame 
remains completely elastic i.e. elastic critical loads have been derived, 
In most cases, an approximate method of analysis has been used, whereby the 
actual loading on the structure is replaced by one consisting of forces 
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applied at the joints only. As long as no lateral loads are present, the 
members of the frame are then subjected to axial load only before buckling. 
This approximate method has the distinct advantage of reducing the amount of 
Calculations required to determine the critical load, and is sufficiently 
accurate for many but not all cases. More recent investigators have used 
exact methods of analysis, whereby the effects of primary bending moments 


and associated deformations are included, 


A recent bulletin(4) published by the Welding Research Council con- 
tains a tabular summary of known solutions for elastic critical loads for 
frames subjected to axial forces only (approximate method) and for frames 
subjected to axial forces and primary bending moment (exact method). In 
addition, a very extensive bibliography on the subject of frame stability 
is presented. Some of the more important of these solutions plus other in- 
vestigations of interest are briefly discussed in the following two sections, 
(A) Frames Carrying Axial Forces Only 

Solutions for the elastic critical loads of frames using the approx- 
imate method of analysis have generally utilized some modified form of the 
slope-deflection equations or moment-distribution procedure, although sol- 
utions based on energy principles have been obtained, Since all members 
are subjected to axial force only, uniform yielding may be assumed for every 
section along each member; hence, an estimate of the inelastic buckling load 
can be obtained from these solutions by applying a tangent modulus modificat- 
don. Lu(5) presents an excellent summary and brief discussion of these 


methods, 


Bleich(6) and Timoshenko(7) include solutions for the elastic critical 
load of simple frames with either fixed or hinged column bases. Merchant 
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and others(8) (9) (10) obtained similar solutions for more complicated multi- 
story, multibay structures, using type solution methods suggested by Bolten, (11) 
and stability functions developed at the University of Manchester by Livesly 
and Chandler.(12) wooda(2) calculated various "deteriorated critical loads" 
of a multistory frame as part of an investigation of his concept of a deter- 
iorated elastic structure for stability analyses. Goldberg(13) investigated 
the critical loads of complete one-story building frames, where the roof is 
considered as an elastic restraint acting between the frames, This repre- 
sents a more realistic approach to the actual situation. Galambos(14) 
demonstrated that only a small degree of restraint at the base of the col- 
umns increases the critical load considerably above that for hinged-base 
colums. A.S.C.E. Manual Ne. 41(15) contains graphically the solution to 
the elastic critical load of an analogous portal frame that was used to for- 


mulate a design rule for the new A.I.S.C. Code. 


The theory of stability for structures whose members are subjected 
only to axial forces has been fully developed, and numerous methods for 
determining the critical load of various types of frameworks are now avail- 
able. For complex structures, however, the numerical work involved in the 
analysis is often excessive, and simplifications in this field are required 
so that the required computations can be readily performed in a design of- 


fice. 


(B) Frames Carrying Axial Force and Bending Moment 

Comparatively few solutions exist in the literature dealing with the 
elastic stability of frames subjected to ae which produce primary bend- 
ing moments as well as axial forces in the members before buckling. The 


calculation of critical loads for frames using the available solutions is 
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often very tedious, and accurate numerical results have been obtained only 
for single bay, single story frames. Solutions to the same problems in the 
inelastic range are non-existent except for an initial attempt outlined in 
a thesis by tu(5), Empirical approaches still remain as the only method of 
estimating a true elasto-plastic failure load for frames subjected to axial 


force and bending moment. 


Bleich(6) gives a summary of the method followed by Chwalla in 1938, 
for a hinged rectangular frame subjected to two equal and symmetrically pla- 
ced loads on the beam. Chwalla used the classical approach of integrating 
a system of differential equations and satisfying the associated boundary 
conditions. More recently, Masur and others(16) succeeded in extending the 
Slope deflection and moment distribution methods to the analysis of the 
stability of frames carrying primary bending moment. Using similar means, 
tu(5) obtained exact solutions for the elastic and inelastic critical loads 
of a hinged base rectangular frame subjected to a uniform load on the beam 
and two equal concentrated loads over the colums,. This particular frame 


was chosen to simulate the lower story of a multi-story building. 


Preliminary but unpublished work at the University of Alberta in 1963 
has indicated that the existence of horizontal as well as vertical loads on 
a simple portal frame results in substantial sway deflections. This work 

thus suggests the use of a so-called "large deflection theory" to obtain 


more accurate elastic critical loads for this combination of loading. 


The increasing use of electronic digital computers has served to 
further the use of existing methods in this field of endeavor, This thesis 
employs the use of such @ computer in order to obtain elastic critical loads 


for building frames described in the next chapter, 
10 
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CHAPTER II 
SCOPE OF PRESENT INVESTIGATION 


Published results of investigations using approximate or exact 
methods of analysis cover cases where the frame is subjected to vertical 
load only. The investigation contained in this thesis was initiated to 
obtain exact solutions for the elastic critical loads of a simple rectangu- 
lar frame under the action of a combination of vertical and horizontal 
loads. Also considered in this thesis is the effect of using a large deflect- 
ion theory, and the influence of partial base fixity on the elastic critical 


loads, 


More specifically, the following frames, loadings, and ranges of 
variables were investigated: 


(i) Frame with Pinned Base 


wle= (Fe). 


6). hii RP = B fox 


N= 1,/L, 
REI = EI, /EI, 


FIGURE 2-1 


Small deflection and large deflection theories were investigated 
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for this frame and loading throughout the following ranges of variables 


(S.D. - small deflection, L.D. - large deflection) 


For all combinations above, RP values of 0, 0.01, 0.04, 0.07, and 0.10 


were used. 


(ii) Frame With Fixed Base 


wl_= XP.) 


@F)p ——= 
RP = B/ x 


N= L/L, 


REI = EI, /EI, 


FIGURE 2-2 


Small deflection and large deflection theories were investigated for 


this frame and loading throughout the same ranges of variables as for the 
hinged-base frame. 


(iii) Frame With Partial Base Fixity 


The effect of partial restraint to rotation of the column bases on 


the elastic critical loads was simulated by inserting a restraining beam be- 


tween the bases of the colums, as follows: 
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a RP = B/a 


N= L/L, 


REI = EI, /EI, 


FIGURE 2-3 "= EL,, /El, 


The small deflection theory only was investigated throughout the 


following ranges of variables: 


For all combinations above, i values used to plot graphs were n= 


0.0, 0.20, 0.50, 1.0, 2.0, 4.0, and 8.0. 


As can be seen, the large deflection theory for parts (i) and (ii) 
and the whole of part (iii) are investigated only at the extremes of the com- 
plete range of variables outlined in part (i). This approach was followed 
since an indication only of these results was considered necessary for pract- 


tical use, 
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CHAPTER III 
METHOD OF SOLUTION 


In order to obtain exact solutions to the situations outlined, the 
effect of primary bending moments and associated deformations must be in- 
cluded. Since it is not possible with the techniques available to obtain 
unique mathematical solutions to problems of this kind, an iterative sol- 
ution using the slope deflection equations modified to include the effect 
of axial load was chosen. The use of the University of Alberta's I.B.M. 
1620 electronic digital computer made possible the use of a method such as 


this. 


The necessary slope-deflection relationships may be developed by 
referring to Figure 3-1, in which a typical frame member connecting joints 


A and B is shown. 


W 5 original position 


CRE, GU 
eae. A ali alae all 


FIGURE 3-1 


The typical member A-B is subjected to an axial force P, which is 
considered positive in compression, end moments M,, and M,,, end shears, 


and lateral load W between the ends. The joint rotations are 0, and 6,, 
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and the bar rotation is @, when the member is in the deflected position. 
End moments and rotations are defined as positive clockwise. The bending 
moments and rotations are defined as positive clockwise, The bending moments 


and joint rotations are related by the slope-deflection equations as follows: 


My, = Ky, (0, + C,,0,) - K,, (1 + G,,) Tan @ * FEM 

May ® Ky (0, + C,,9,) - K,. (1 + C,,) Tan 0 + FEM 
In general, M,, = K..(0, +00.) = K..(1+C,,) Tan 9 + FEM where N de- 
notes near end and F denotes far end, For relatively small rotations, Tan 
@ effectively becomes 0, In the slope deflection equations, K is the stiff- 
ness factor, C is the carry-over factor, and FEM is the fixed-end moment for 
the lateral load on the member. If the modulus of elasticity E is constant 


and the member ie prismatic, then K,, = K,, = K # SEI 


The stiffness coefficient S and the carry-over factor C are functions 
of the axial load P in the member, as shown in Figures 3-2 and 3-3. These 
functions have been tabulated(12) for a wide range of axial loads, and are 
briefly developed in APPENDIX A. For the case of no axial load, the S and 


C functions reduce to the conventional values of 4.0 and 0.5 respectively. 


The fixed end moment is a function of the magnitude of the axial 
load P and the magnitude and distribution of the lateral load W. For a 
prismatic member, 

FEM,, = -K(Oy, + CO,, ) 
where 0,, and 0, are the rotations at the ends of the frame member when 
considered as a simply supported beam-column subjected to axial and lateral 


load. For a member subjected to a uniform lateral load w, and an axial 


15 


pene atbaed siT ‘Yeaiwloolo evivleod 64 bontteb ha enordaaer fin 


iewekict es enolistipe aolics [teb-o7q ak tit vd bejefet es onokintot 3 


a —" 


a 
Mit * © ant ( 9? + £) oa - oe oa” * ;®) ad? «na ad 


‘oa 
oa 
~ 


i) 
2 


Wi = 6 ast ( uw +5) of © aa? + 59) itl ad ; 


+ . . " / ' AA : 
eb A erect wi = & rat i me s rf) A ~ \.9 ey” * a) = = hts 


a 


ay 
a ' , an 
na? ,atolbidnicr ifm vievitsle? to% bree tet estoneh T bam fem 


‘ : 
oe 


= 


~tiise eff et 1 ,eaoltsyvpe motieaiteh egofe si? al .9 nennoed: gt 
“Ot Snenem boe-~textt af? ef SAT bas ,wotorlt tevoeyrtao edt db O om 

tnadvenco «ft 7 yelolieale Yo svlobam erie TI = .tedmee ed? Ao hood Ly 
(SA A» yt opdt ,oldaceing ot aedae 


i) wr ag a *) a 


smolfont'! o1g D todos? sevo-ritso ect baa 2 tnetortteoo agent 


seent .£-@ bao S~E exurti nt mrofs ea ,sodmem edd mt 4) Been 
vy) a aT 


a ‘ " = sre r . Le Py 
ets bas .sbaoi Isixs to syiusi shiw g tot \>- ted afudad need evad - 
& es 
ti. eun 
i 
bas 2 ont ,beot Lakie on to sess edd 207. .A LIGMESGA ae 


«Ylevtiosqse: ¢.0 bag 0,4 to eoulsy fancisaevace eff oF 1 


lefon ett be sbhudtngsn end to aofgomt 5s et 


6101 4* deol fetsdal ent Yo aoksedinteath bas Obs 


feu todmen enert ds to abne edd Ja enolListor si 
Lapetal baw Iaixe ot be toefdwe. fiuss.Loo-rrasd et ot a 


faixe os fra, gw heat lanodal meat. Bits & od 


compressive force P, the expressions for the fixed end moment have been deve-~ 
loped in APPENDIX B and are shown graphically in Figure 3-), For the case 
of no axial load, the fixed end moment reduces to the conventional value of 


0,083 wL, 


The use of the slope deflection equations in the structural analysis 
of a building frame is illustrated fully in APPENDIX C, The slope-deflec- 
tion equations are written for each joint of the frame, with the joint 
rotatiors and sway rotations as unknowns, By writing the equilibrium equa- 
tion 2M = 0 at hinged joints, it is possible to express these joint ro- 
tations in terms of others, thus reducing the number of unknowns, Writing 
the equilibrium equations <M = 0 at the remaining joints and ZF = 0 for 
the frame as a whole, the governing equations for the solution of the joint 
and sway rotations of the frame are obtained, The equations will contain 
a load parameter of the frame which must be given a value in order to solve 
for the rotations, Conversely, if one of the rotations is specified in 
value, the load parameter and the remaining rotations can be solved, This 
approach is used in this thesis to solve for the elastic critical loads of 


the frames outlined in CHAPTER II, 
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CHAPTER IV 
PRESENTATION AND DISCUSSION OF RESULTS 


4.1 GENERAL 

In present-day American design practice, it is quite common to 
assume column bases of frames as hinged, However, it has been shown 1h ) 
that presently used column bases offer a considerable restraint to base 
rotation, This restraint is enough to raise the critical load of some frames 
considerably above that of a hinged=base frame, Furthermore, the effect 
of cladding, composite action, infilling, and bracing in actual buildings 
is not considered in this or most other stability analyses, but is def- 


initely beneficial towards preventing frame instability. 


The results of this investigation are presented graphically on the 
following pages, For clarity, the results are divided into three separate 
sections, Sections h-2 and }}-3 contain the results for the hinged-base 
frame and fixed-base frame respectively, Section h-l contains the results 
for the frame with partial base fixity, this case representing various 


degrees of base fixity between the extremes of hinged and fixed bases. 


<2 FRAME WITH HINGED BASES 
The relationships between the load parameter < and sway rotation 
@ are shown in Figures -3 to )-26 for the five values of horizontal load 


considered, The results of the "large deflection theory" are given following 


those of the small deflection theory. 
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For the case of no horizontal load, the maximum vertical load the 
frame can support if it remains elastic is 0.5 brad ne. which occurs with an 
infinitely stiff beam member, This can be deduced by considering one of 


the column members in its buckled configuration under this condition, as 


in Figure -2, 


— 
\ 
ae | 
: | | rotation fixed 
> | translation free 
oi FIGURE h-2 
mi rotation free 


translation fixed 


The effective length of this column member is 2L, making the elastic crit- 


ical load equal to 


WEI = 0,25 TEI = 0,25 P. 
(20) pa 


The frame has two such column supports, and hence its elastic critical 


male 
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load becomes 2 x 0,25(P, ), = 0,50(P, ), 


It can be seen by considering the results of the small-deflection 
theory that the maximum vertical load sustained by this frame with no hori- 


zontal load present does in fact approach Oyo ). as the beam stiffness 


c 
increases, Conversely, when the beam stiffness is relatively small, the 
critical load drops considerably since the effective lengths of the colums 


are free to increase above 2L, 


The existence of even a small horizontal load causes a significant 
increase in sway rotations, with larger sway rotations being evident for 
larger horizontal loads, The effect of horizontal loads on critical 
vertical loads was difficult to assess quantitatively, since a definite 
critical load was not always reached within a range of sway rotations less 
than 1.0, If it is assumed that the results are reliable throughout the 
range, and if the critical load is defined as the maximum load occurring 
in the range of 6 = 1.0, then it may be stated in general that the exist- 
ence of a horizontal load does not decrease the critical load substan- 


tially, but has more effect on sway rotations, 


The results of the "large-deflection theory" for only the extremes 
of the variables presented for the small-deflection theory are given in 
Figures -19 to h-26, For all cases the total vertical load at a given 
sway rotation as predicted by the "large-deflection theory" is larger than 
that predicted by the small-deflection theory, This phenomenon is most 
pronounced at large sway rotations, while for very small sway rotations, 


the results are identical for both theories, 


Tt should be noted that both the large and small deflection theories 


are based on the well-known approximate relationship M = EI d‘y/dx*, which 
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describes the equilibrium configuration of a member, considering flexural 
deformations only, However, this equation is valid only if the difference 


between the total length along a member and its corresponding chord length 


is negligible, The exact value for curvature is d°y/dx* from which 
1 + (dy/dx)*}” 


2 [ 

M = EI d-y/dx* ,2 if an error of 5% is allowed in this equation, then 
1 + (dy/dx)*|% 

it can be shown that dy/dx must not exceed 10,)° or 0,1815 radians, An 


investigation of the data results from the computer programs has indicated 
that an error of this magnitude will occur for sway rotations of the frame 
exceeding © ~ 0,20 radians = Tisbs Therefore, the results of both the 
small and "large" deflection theories mist be considered inaccurate for 
sway rotations beyond 9 = 0,20 radians, 

The above discussion indicates that a true large-deflection theory 


should be based on the exact relationship M = EL d°y/dx* The large- 


1 + (dy/dx)'|% 
deflection theory used in this thesis must therefore be called a "pseudo 
large-deflection theory", since it considers only part of the difference 
between small and large deflections, as discussed in APPENDIX C, 

It has been decided to treat the results of the small deflection theory 
as a reasonably conservative estimate of the elastic critical loads of this 
frame, and to use them in forming the following general observations: 

(1) The elastic critical load decreases with increasing horizontal 
load, all other variables held constant, The percentage reduc- 
tion from the case of no horizontal load is about the same for 


all combinations of \ and REI investigated, 


(2) The elastic critical load decreases with increasing values of \, 
the variables RP and REI being held constant, 


(3) The elastic critical load increases with increasing values of REI, 
the variables RP and \ being held constant, i.e., the critical 
load increases with increasing beam stiffness, 


A summary of the elastic critical loads obtained is given in tabular 


form following the graphs, 
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The relationships between the load parameter « and sway rotation 
@ are shown in FIGURES h-29 to h-52 for five values of horizontal load 
considered, The results of the large-deflection theory are given following 


those of the small deflection theory. 


For the case of no horizontal load, the maximum vertical load the 
frame can support if it remains elastic is 2,0(P, )., which occurs with an 
infinitely stiff beam member, This can be deduced by considering one of 
the column members in its buckled configuration under this condition, as 
in FIGURE )-28, 


rotation fixed 
translation free 


FIGURE 28 


Fer LENGTH =L 


rotation fixed 
translation fixed 


The effective length of this column member is 1,0 L, making the 
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elastic critical load equal to 


TEI = 1,0 P, 
L? 


The frame has two such column supports, and hence its elastic critical 
load becomes 2x 1,0 (P,). = 2.0 (P.).. 

It can be seen by considering the results of the small deflection 
theory that the maximum vertical load sustained by this frame with no hori- 
zontal load present does in fact approach 2,0 (P,). as the beam stiffness 
increases, Conversely, when the beam stiffness is relatively small, the 
critical load drops considerably since the effective lengths of the columns 


are free to increase above L, 


Difficulty was encountered with some of the cases in this series, 
For some of the cases where ) was relatively large and REI relatively small 
(A= 3,0, REI = 3,0 and 6,0, and A= ,0, REI = 4,0, 8.0, 12.0), the com- 
puter results were very erratic as the critical load was approached, In 
some cases, it was possible only to complete the lower portions of the 
load deflection plots, For the remaining graphs, the results were extra- 
polated from the partial plots and the complete results of other graphs, 
Any extrapolated curves are shown as broken lines instead of solid ones, 
The reasons for this difficulty are not definitely known, but it is believed 
the mathematical solution becomes unstable as the critical load is approach- 
ed for these particular combinations of variables, It is assuring that these 
difficulties have occurred in the less practical portion of the range of 
variables considered, since it is probable that larger beam stiffnesses 
would be used when ) is large, 

As in the case of the pinned base frame, the existence of even a 
small horizontal load causes a significant increase in sway rotations, with 
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larger sway rotations being evident for larger horizontal loads, Again, 

the effect of horizontal loads on critical vertical loads was difficult to 
assess quantitatively, since a definite critical load was not always reached 
within a range of sway rotations less than 1,0, On the basis of the results 
of the small deflection theory, it may be stated in general that the exist- 
ence of a horizontal load does not decrease the critical load substantially, 


but has more effect on sway rotations. 


The results of a "largeedeflection theory" for only the extremes of 
the variables presented for the small-deflection theory are given in Figures 
yelS to h-52, As for the pinnedebase frame, the large-deflection theory’ 
produced higher loads than the smalledeflection theory in the range of large 
sway rotations, whereas the results are identical for very small sway rota- 


tions, 


As discussed for the hinged-base frame, the use of the equation 
M= EId“y/dax’? introduces errors to both the small and large deflection 
theories used in this thesis, Allowing an error of 5% between this equation 
and the corresponding exact equation, and investigating the computer results 
for this frame, it was found that an error of this magnitude will occur in 
some cases for sway rotations as small as 9 = 0,05 radians, Therefore, 


the results of both the large and small deflection theories must be con- 


sidered inaccurate for larger sway rotations, 


Considering the results of the small-deflection theory as a reason- 
ably conservative estimate of the elastic critical loads of this frame, the 


same general observations as for the pinned base frame can be made: 
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(1) The elastic critical load decreases with increasing horizontal 
load, all other variables held constant, The percentage reduction 
from the case of no horizontal load is about the same for all com- 
binations of A and REI investigated. 


(2) The elastic critical load decreases with increasing values of ) ’ 
the variables RP and REI being held constant, 


(3) The elastic critical load increases with increasing values of REI, 
the variables RP and ) being held constant, i.e., the critical 
load increases with increasing beam stiffness, 


A summary of elastic critical loads obtained is given in tabular 


form following the graphs, 
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hel) FRAME WITH PARTIAL BASE FIXITY 


RP = &/. 
= L/L, 
REI = EI,/EL, 


" = (EI/L),,/(EI/L), = EL,,/EI, 


FIGURE he53 


One more variable is introduced for the frame with partial base fixity 
compared with the previous two cases, The relationships between the elastic 
critical loads and the restraint parameter » were plotted and are presented 
in Figures 5), to )-61 for the two values of horizontal load considered, 

In plotting these graphs, the critical load used was the greatest value of 
the load parameter « obtained for a range of sway rotations less than or 
equal to 1,0, 

When y = 0, the elastic critical load is the same as for the hinged 
base frame, As 1 increases above zero, the elastic critical load increases 
rapidly at first, but levels off to approach asymptotically the value of 
the elastic critical load for the corresponding fixed-base frame as 1 ap- 
proaches infinity. 

It can be observed for all combinations of variables investigated 
that only a small amount of base restraint increases the critical load 
considerably, This behavior is most pronounced when the beam stiffness 


is large. If the base restraint is equal to the stiffness of the top beam, 
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h-5 DESIGN EXAMPLE 

To illustrate the use in actual design of the elastic critical loads 
derived in this thesis, the following design example is given, Using 
Merchant's empirical approach (as discussed on page 7), an estimate of the 
true failure load of this example frame will be obtained. Preliminary analysis 


indicates the horizontal loadiwill be about 10% of’ the total vertical load, 


| ee ti EB = 30x10" ksi 
cn Sos %y = 33 ksi 
ca - Neglect frame weight, 
| 30' | 
a 
FIGURE )-62 


(A) Ultimate Strength - Simple Plastic Theory 


=) 3 7 
uw va Sees 
2 ‘ “A 28 
Mech,#1 Mech, #3 
3 w(150) = 2M, 0 30 w(7,50) = 6M, © 3 w(150) + 30w(7.50) 
w= O,O4L5 My, we 0,0266 Mp = 6Mp, 0 


w= 0,0222 Mp 


Y% 
Mech, #3 governs. Hence, ultimate load w= 0,0222(83) = 1,8), 
(B) Elastic Critical Load 
\ = 30/15 = 2,0 REI = 273/110 = 2.5 RP = 0,10 


interpolating between figures +7 and 4-8, “erie 0,345 
e/¢ 


w= 0.345 (Pp = 0.35 I 30200, \.110) = 11.53% = Ure 
30( 1 
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(C) Actual Load-carrying Capacity 
The best available estimate of the true failure load can be ob- 


tained by the application of Merchant's empirical equation (discussed on 


page 7) * 


where wu, = true failure load 
Ug = elastic critical load 
wy = Simple plastic collapse load 


Therefore} oneldetieltinn 1 
Wie Lusas Toy 


We = 1.59% = true failure load. 
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CHAPTER V 
CONCLUSIONS AND RECOMMENDATIONS 


The following general conclusions regarding the results of this 
investigation can be made: 

(1) A slope-deflection analysis using an electronic digital com- 
puter is a suitable method for solving stability problems for frames sub- 
jected to axial loads and primary bending moments, However, the mathemat- 
ical solution became unstable for the frame with fixed bases when the length 
of beam was large relative to the column length, and when the beam rigid- 
ity was small relative to the column rigidity. Similar difficulties may 
be encountered in future stability investigations of other building frames, 

(2) A fixed-base frame of the type considered in this thesis 
will sustain a elastic critical load about four times greater than the 
corresponding hinged-base frame, for a given value of horizontal load, 

The load-deflection curves for these two frames are similar except for 
this ratio of loads, When horizontal load is present, both frames ex- 
hibit large sway deflections before the critical load is reached. 

(3) Conclusions (4), (5) and (6) are applicable to both the hinged 
and fixed-base frames. 

() The elastic critical load decreases with increasing horizontal 
load, the frame dimensions and ratio of member rigidities remaining con- 
stant, The percentage reduction from the case of no horizontal load is 


about the same for all combinations of variables investigated, 
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(5) The elastic critical load decreases at an increasing rate as 
the length of beam increases with respect to the length of columns, the 
horizontal load and ratio of member rigidities remaining constant. 

(6) The elastic critical load increases at a decreasing rate as the 
beam rigidity increases with respect to the column rigidity, the horizontal 
load and frame dimensions remaining constant. 

(7) The use of a so-called "large-deflection theory" results in 
higher load parameter values being obtained for a given sway rotation than 
those obtained by the small-deflection theory, This result is most pro- 
nounced for large sway rotations. 

(8) The presence of even a small amount of base restraint increases 
the critical load considerably above that of a hinged base frame, This 
amount of base restraint is available under presently used column bases 


that are commonly assumed in design to be hinged bases, 


The following fields of interest are suggested for future work in 
the elastic stability of frames subjected to axial loads and primary bend- 
ing moments: 

(1) The critical loads of portal frames under the action of uniform 
vertical and horizontal loadings may be investigated. 

(2) The critical loads of multi-bay, multi-story rectangular build- 
ing frames subjected to horizontal and vertical loads may be investigated. 

(3) The critical loads of gable frames subjected to various com- 
binations of loading may 53 investigated. 

(4) The effect of gusset plates in the corners of these frames may 
also be investigated, Stability functions for members with gusset plates 


2 
at the ends have been developed ) 
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APPENDIX A 


STIFFNESS AND CARRY-OVER FACTORS 


Consider a member of length L and flexural rigidity EI. 


Define forces P and F, moments M, and M,, and end rotation 0. 


positive as shown in figure A-1l. 


F = Ma+ Me 
Scalar eee 


| 
) SS Oct ace 2 FIGURE A-1 


By definition, the stiffness factor K = SEI/L is the moment 
required at the near end of a member to rotate the near end through 
aunit angle. The carry-over factor C is the ratio of the moment 
induced at the far end to the moment applied at the near end. The 
stiffness coefficient S and carry-over factor C are then given by 
the relationships M = SEIO,/L and M, = CM,, and will be calculated 


for a member with one end fixed, as shown in Figure A-l. 


The differential equation expressing the curvature of this 


member is based on the well-known relationship EI d’y = M and hence 
ax? 


becomes: 
EI d’y = +M, ~ (My + Ma)x - Py 
dx L 


The general solution to this equation is: 


A-1 
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y Te ara M, - (My + Ma) x) where k = 
P L 
A and B are constants to be determined by the boundary conditions. 


Applying the boundary conditions of y= 0 at x 


the general solution becomes: 


ay See Conx))M, + |x - Sin ler) M ; 
ar Sin kL L Sin kL 
To obtain the carry-over factor C, it is necessary to obtain 


the value of M,/M, for the slope at A equal to zero, as follows 


kL Cos k(L-x) + 1] M, s+ [2 ~ Ke Cos tex] Ma 
x = all Sin kL Sin kL 


dy =O meri Sefee te +1] + [2 - 
dx/x = 0 Lk? EI 


Ms} 

Sin kL Sin kL 

oo |C=M = - Sin where % = kL 
Me Sin @ - G Cos @ 


To obtain the stiffness factor K, it is necessary to find the 


Slope at B and set this equal to unity while the slope at A remains 
zero, as follows: 


When dy = 0 at x = O, the general formula for slope becomes 
dx 


dy =__1 _([-_kL Cos_k(L-x) + 1 a - Sin kL |e 
Lk* EI Sin kL Sin kL - kL Cos kL 


+ - - kL Cos id Me 
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The above € and S functions were derived for the case of axial 
compression. Similar derivations for the case of axial tension yield the 


following expressions for C and S: 


If the axial load is zero, it would be expected that the C and 
S functions would reduce to the conventional values of 0.5 and ).0 re- 
spectively. However, the substitution of g = O causes all of the above 
expressions to reduce to zero. This apparent discrepancy can be overcome 


by the application of L'Hospitals rule, which reads as follows: 


Lim f(x) = Lim tx) 
Bayip (x 490 oi (x 


Therefore, successive differentiation of the expressions for stiff- 
ness coefficient and carry-over factor with respect to g and the substit- 
ution of g = O will result in C = 0,5 and S = 4.0, which are the correct 
values when the axial load is zero. 

It must be noted that the stiffness and carry-over factors are 
based on the equation M/EI = d*y/ax*, This equation is valid only when 
the difference between the total length along a member and its correspond- 
ing chord length is negligible. Therefore, these stiffness and carry-over 


factors are not strictly correct for a true large deflection theory, 
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APPENDIX B 
SLOPE-DEFLECTION EQUATIONS 


The resultant end moment on a frame member under the applied 
loading is made up of the following components: 
(A) Moments due to rotations of the ends of the frame member 


(B) Moments due to lateral displacement of the ends of the 
member 


(C) Moments due to effects of lateral loads applied to the 
frame member under consideration 

The magnitude of the axial load present in the member affects 
the resultant end moment and must be considered in any atability ana- 
lysis. It is a well-known fact that the effects of axial load and 
moment cannot be evaluated separately and then combined to obtain a 
correct answer, The true resultant end moment must be obtained by 
considering the effect of axial load on each of the component moments 
above, and then superposing these components by direct addition. 


This is the basic principle of the slope-deflection equations, 


The component moments making up the resultant end moment will 
now be considered separately in the following sections. Note that 
the slope-deflection equations neglect the effects of shear or axial 


deformations, which are relatively small compared to flexural deforn- 


ations, 


(A) Moments due to end rotations 


Consider a member of length L and flexural rigidty EI. Define 
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forces P and F, moments M, and M,, and end rotations 0, and 0, pos- 


itive as shown in figure B-l,. 


F- TatMe 
FIGURE B-1 


As given in APPENDIX A, the general solution to the differ- 


ential equation defining the equilibrium configuration of this member 


is 70 211 Eas Mex) | ¥, *|z - Sint i | Me appar en 2 


Sin kL L. Sin: kL 


The general expression for slope becomes 
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Cross-multiplying, 


i 


Q,LK°EI Sin kL = (Sin kL - kL Cos kL)M, + (Sin kL - kL)M, @ 


(Sin kL - kL)M, + (Sin kL - kL Cos kL)M, @) 


Q,Lk*EI Sin kL 


From} 


M, = 9 Lk EI Sin kL = (Sin kL - kL)M 
Sin kL - kL Cos kL 


Substituting into (2)and manipulating terms, 


M, = KL[__Sin Ml. — ML Gos Wl __)EI 0 
2 —- 2 Cos kL - kL Sin kKLI|L 


+ kL[__ Sin kL - kL Cos kL kL - Sin kL |B 
2 - 2 cos kL - kL Sin kL||Sin kL - kL Cos kLIL 


or | M, = KO, + KCO, = K(0, + CO) = BEG + C0, ) 
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This expression was derived for the case of axial compress- 
ion. A similar derivation for the case of axial tension will yield 
the same expression if the corresponding carry-over and stiffness 


factors for this case are utilized. 


It should be noted that the axial load P as defined is the 
load in the member lying on a line joining the ends of the member. 
For a frame member subjected to sway, it is often more convenient 
to calculate the axial load as the load on a line joining the ends 
of the member in the undeflected configuration, For small sway def- 
lections, these axial loads are nearly equal, and so little error is 
introduced. However, for large deflections, the axial load lying 
along a line joining the ends of the member in the deflected con- 
figuration must be used. This concept is illustrated in the analysis 


of the portal frame in APPENDIX C. 


(B) Moments due to lateral displacement 


Consider the same member subjected to a sway deflection , as 


shown in figure B-2, 


FIGURE B-2 
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The differential equation defining the equilibrium configur- 


ation of this member reads as follows: 


The general solution to this equation is: 


P oo heme 


arlama mate = ase, & [2] where k =|[P_ and A 
L EI 


and B are constants to be determined by the boundary conditions. 
Applying the boundary conditions of Y = 0 at x = 0 andx=L, 


the general solution becomes: 


= M pin tot + Sin lor - Cos kor - 2x + 1] 
k*EI|Tan kL Sin kL 


ci 


The general expression for slope them becomes: 


dy = _M Petes tes + Cos tox +k Sin lox ~ 2 


dx k*EI| Tan kL Sin kL L 


At x = 0 dy = -Tan 0 = M | cee fiw 2 2] 


dx k*EI|Tan kL Sin kL L 


Cross-multiplying and manipulating terms, 
M=@/_ Sin @ - G Cos |Ei> +524 lane 
2-2 Cos @ - G Sin ZI|L Sin % - g Cos 
where = kL ,",|M = K (14C) Tan 0 = SEI (14C) Tan 0 
L 


This expression is derived for the case of axial compression. 
A similar derivation for the case of axial tension will yield the 
same expression if the corresponding carry-over and stiffness factors 


for this case are utilized. 


It should be noted that the axial load P as defined is the 


load in the member lying on a line joining the ends of the member, 
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For small sway deflections, little error is introduced if the axial 
load used is that lying on a line joining the ends of the member 
in the undeflected configuration. Furthermore, for small sway def- 
lections, Tan 9 is effectively equal to 9, and the moment due to 


sway can be taken as M = K(14C)O, 


(C) Moments due to lateral load 

This component moment is simply the fixed-end moment for the 
member, Since end rotations are considered under part (A). The fixed 
end moment is a function of the magnitude of the axial load and the 
magnitude and distribution of the lateral load. 

In general, for a prismatic member, 

FEM, ,= -K(0,, + CO, ) 

where 0,. and 0,, are the rotations at the ends of the frame member 
when considered as a simply-supported bean-column subjected to axial 
and lateral loadings. This relation is an application of the results 
of part (A), since the ends of the member are rotated from their sim- 
ply-supported position to that of no end rotation in order to get the 


end moments, 


It is apparent that expressions for fixed end moment depend 
on the type of lateral load distribution. Since uniform lateral loads 
only are considered in this thesis, this case only will be considered 


here, 


For a member subjected to uniform lateral load w and no axial 


load, as assumed in conventional frame analysis, 
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Consider further a member subjected to uniform lateral load 


uy and axial compressive force P, as shown in figure B-3, 


A 
p —_ Carter, FIGURE B-3 
7 (OEE Eero pe 
‘ | As bs 
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The expressions for @,, and @,, can be derived by considerat- 


ion of the differential equation (7) 


oat aes 


The general solution to this equation is: 


Y = A Sin kx + B Cos kx + Cx + D + wx’ =where k = |\/P_ 
2P EI 


and A, B, C, and D are constants to be determined by the boundary 
conditions, 


dy = Ak Cos kx - Bk Sin kx + C + wx 


ax P 
dly = -Ak’ Sin kx - Bk’ Cos kx + w 
dx P 


Applying the boundary conditions of y = O and d°y/dx*= O at 


both x = 0 and x = L, the general solution becomes 
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or @, “wE/2 Tan g =i w= 6. where 9 = kL 
2P | 9 
Therefore, for a member subjected to a wiiform lateral load 
wand an axial compressive force P, 


FEM,, = -SEL [WL a-0)/2 Tan § ~ 1 
L |2P Gi we 


or |FEM,, = -S(1-C) [2 Tan g ~ lwi? 
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A similar analysis for the case of a member subjected to a 
uniform lateral load w and an axial tension force P will yield the 
following expression for fixed end moment, using the corresponding 


stiffness and carry-over factors for this case: 
FIM. = = Sa - 4 Tanh g) wl? 
Ne la 5 


These expressions are strictly correct only when the axial 
load used is that lying on a line joining the ends of the member 
in the deflected configuration, the lateral load used is that perpen- 
dicular to this stated axial force, and the rigidity EI is uniform 
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APPENDIX C 
STRUCTURAL ANALYSES OF FRAME 


Consider the portal frame of figure C-1 with the hinged colum 
bases, and assume that the deformations in the equilibrium configur- 


ation of figure C-lare small but finite. 


| urLy= off C 


FIGURE Cel 


The slope deflection equations are written as follows: 
Mig = Kyy (0, + Cyp OQ) ~ Kyg (1 - Org) © 

Ma = Ky (9, + 0,9) -K,,(1 + C,,) © 

=K,. (0, + 0,,0,) - FEM 

Meg = Ky, (OQ, + Coy 0,) + FEM 

Men = Kop (OQ + Gey Op) — Key (1 + Coy ) 0 

Mi = Ske G + Cop 9% ) - K,, (1 + C,,) @ 

Summing moments about A, 


Kyn (0, + Crp Og) = Ky (1+ Cy) O or O = (14+ 0) O - C0, 
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Summing moments about B, 


Kyp (On + Cp) - Kg (1 + CO, 0 + Ky, (0, + 0,0.) - FEM 


it 
© 


Substituting ©, from above and collecting like terms, 


- 2 
[Eup (1 - Gre ) + Ky] O + [K,. Cy, | = [Kip (1 - fy )]@ = FEM 
Summing moments about D, 
Key (0, + Cy O,) = Ky (14Cep JO or 0, = (14g EO - CO, 


Summing moments about C, 
Kee (OQ + C,, Og) + FEM + K,, (0, + C0.) -K,, (1 +O, ) @ = 0 


Substituting 9, from above and collecting terms, 


[Kee Cre] G5 + [Ky (1-02, ) +k] - [K, (1-0, Joe -Fm © 


The equilibrium equations for beam BC are: 


+) 2M= 0 = My, + Mop + of(P, ), be : 


V, = (Pele + (Mpc + 
2 Ls 
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V, = <6 Ye = (Mee Mes) 
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FIGURE C-2 


The equilibrium of column AB is as follows: 


+)2Mj= 0 = My+ WA + BL, 


H = -V,4 - Mpa 
A E; 


FIGURE C-3 
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The equilibrium of colum CD is similar and results in: 


H, = - WA - My 
Ly 


Summing the horizontal forces for the whole structure, 

H, + H, = @(B ), 

YA +My + VO +M, = - O(R Jae 

Substituting for V, and V,, 

(Pe), A + May + Mop = - 6(Pe ) Ly 

Substituting for M,, and M,,, and subsequently for 0, and 9, 
and collecting terms, the third governing equation is found to be: 

o(Pe A+ [Kg (1 - Cig )]O, + [K,(2 - C2, Jo, - [K, (2 - cy ) 

+ Kyp (1 - Coy de = -AF, Lb, ® 

Rewrite equations 1} @)and G) substituting: 

K=SEI/L A=L,0 6=EPa P, = TEI/? d=1,/L, 


and FEM = S(1 -C)F wi = S(1-C)F«(P.),L, = SQ-0)Recnt ED 
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transpose terms to get the required form of the governing equations: 


[Sov (1G, ) + Sec|%s + [See ae] % - Su, Coe ) oat Fee ‘ [Ssa¥ (O-Crp ie O 


[Sec Cyc | O, + [Sco ¥ (1-02 ) + Sse] & + [Sec (AnGy, )art Fac] = [See¥ (1-c2, ie @ 
Se 1-04 210 + [5,, (1-02, Jo, + irae [5,02 ) +8,,0-08 Jo @ 
For the calculation of the S, C, and F functions, the axial 
forces in the members must be known, 
The axial force in column AB is effectively: 
Ph, =V, = Ske Mee Moa ) which after substitution be- 


comes: 


Then ?,¢ X= Sac (14Ca¢ )(On+0¢ ) 
can 2 1? Yd 


The axial force in colum CD is: 


Then Pop 7 5 + Sax UsGoe )\Gas0e) 


The axial force in beam BC is: 


Ps, =H, = -YpA = Mcp and by making the appropriate sub- 
C 
stitutions, , = Pec = -(x¥A+ Soc (1#Cec)(On+Oc) |9 - Sep¥A(1-Co, _)(O, ~0) 
(Pe J 2 m* tT? 


The computer program as described in APPENDIX D solves the 
governing equations to obtain a value of « for given values of 9. 
An iterative solution must be used, since the governing equations 


contain the S, C, and F functions which depend on the axial load. 
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The above derivations have assumed small deformations of the 
frame, A so-called large deflection theory can also be developed by 
reference to figure C-4 which shows the frame in its equilibrium con- 


figuration under large deflections, 


FIGURE C-4 


The difference between the large deflection and the small deflec- 
tion theory considered in this thesis is basically two-fold: 

(i) When writing the equilibrium equations for members with 
substantial sway deflections (such as column AB or CD), a reduced 
length LCos® must be used when summing moments, For small deflect- 
ions, LCos@ becomes effectively equal L, and the equilibrium equat- 
ions become that of the small deflection theory. 

(ii) When calculating the axial load for members with sub- 
stantial sway deflections, the load on a line joining the ends of the 
member in the deflected configuration must be used. For small deflect- 
ions, the load lying on a line joining the ends of the member in the 


undeflected configuration can be used with very little error. 


Following this reasoning, the slope-deflection equations for 


this frame are the same as for the small deflection theory except 0 
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is replaced by TanO, as discussed under APPENDIX B, Hence, summing 


moments about points A,B,C, and D results in the same equations as 


the small deflection theory except 9 is replaced by Tan 9. The equil- 


ibrium equations for beam BC remain exactly the same. The equilibrium 


of column AB changes, however, and is as follows: 


+) 2M,= 0 = My + VA + HL.Cosd 


H,= -VaA - Mes 
L,Cos © 


FIGURE C-5 


The equilibrium of colum CD is similar and results in: 
H, = -VpA - Mep 
L, Cos © 
Summing the horizontal forces for the whole structure follows 
the same routine as before, and results in the third governing equa- 
tion which is: 
(Pe) A + [Kyg(1-Chp 1], + [Key (1-Cp, JO - [Kyg (1-Cyg ) + Key (1-0, )] Tan © 


= -8(R ), L, Cos 0 3) 


Rewriting equations (Y2) and G)by making the same substitutions 
as before except A= L,Sin®, dividing by (EI/L),, and then transpos- 


ing terms, the required form of the governing equations becomes: 
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The axial force in colum AB is derived as follows, consider- 


ing figure C-6, 
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After substituting for M,., Mg and M,,, and simplifying, the 
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FIGURE C-6 
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similar analysis for column CD yields: 
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The small and large deflection theories have now been devel- 
oped for ae frame with the hinged bases. The development of the gov- 
erning equations for the frames with fixed base and partial base fixity 
have been completed and are available at the Department of Civil Ingin- 
eering of the University of Alberta. The analyses parallel those of 


the previous sections, 


A summary of the governing equations for the frames considered 


in this thesis is outlined on the following pages. 
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(A) Small Deflection Theory it i FIGURE C-7 
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Frame with Fixed Bases 


(A) Small Deflection Theory FIGURE c-8 
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Frame with Partial Base Fixity 
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It should be noted that the results for all frames are reported 
in CHAPTER IV using the variable REI = (EI),/(EI),. Knowing ¥Yand), REI 


is obtained from the relation REI = )/¥. 
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APPENDIX D 
COMPUTER PROGRAMS 


The use of the University of Alberta's I.B.M. 1620 Electronic 
Digital Computer greatly facilitated solution of the governing equa- 
tions given in APPENDIX C. Since the S, C, and F functions depend on 
the axial loads in the members which are not known until the govern- 
ing equations are solved, an iterative solution is necessary to ob- 
tain the desired degree of accuracy. Values of S, C, and F of 4.0, 
0.5, and 0.083 respectively for no axial loads present were used to 


obtain the initial solution of the governing equations. 


The flow diagram given in figure D-2 briefly describes the 
general sequence of operations that was introduced to the computer 
for the solution of the critical loads of the portal frames analyzed 
in this thesis. Basically, the computer was instructed to calculate 
successive values of the load parameter «for successive increments 
of the deflection parameter 9, The point on the load-deflection plot 
at which deflection increases with no increase in load defines the 


critical load, as in figure D-l. 


CRITICAL LOAD 


FIGURE D-1 
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Read data card containing 
values of the quantities 


for all members, set 
S=4.0 C=0.5 
16 


Solve original set of 
Simultaneous equations 
to get load « 
(4) 
Calculate axial loads in all members 
and check whether tension or compression 
(5) 
Calculate S, C, F functions 
| (é) 
Solve set of simultaneous 
equations to get new load «, 
|) 
If 'newn"' is within 3% of 
"old" &, print out answer as oj. 
If not, repeat 4, 5, and 6 until 
it does, Then go to &. 
(8) 
If «; 4 %;, , end program. 
If not, increment 9 and go 
to 2. 


End 


FIGURE D-2 Flow Diagram for Computer Program 
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The governing simultaneous equations were solved by using a 
determinant type of solution. 1f the number of equations is larger 
than three, it would be more expedient to use some other solution 
technique, such as a Gaussian elimination procedure. As an alternat- 
ive to including this portion in the main program, a computer library 
program is sometimes available and can be called for by the main pro- 


gram. 


In solving for the load parameter «, the program was instructed 
to repeat itself until the error in this value was less than one half 
of one percent. The use of this particular standard of accuracy was 
by choice nominal, but the results have indicated that this is a rea- 


sonable value since satisfactory load deflection curves were obtained, 


The value of deflection used to increment each successive value 
of 9 had a very direct bearing on the running time of the programs, 
For this reason, a relatively small 9 increment was used initially in 
order to adequately define the rising portion of the load deflection 
curve, while for the flatter portion of the curve, the computer was 
instructed to increase the increment value. The initial values of the 
increment used for various values of RP, the ratio of the lateral to 


vertical load, are as follows: 


RP @ INCREMENT (RADIANS) 
0.00 0.05 
0.01 0.02 
0.04 0.03 
0.07 0.04 
0.10 0.05 
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After four points were calculated, the computer doubled the increment, 
and after twelve points, the increment was made three times the ini- 


tial value. 


A preliminary program was run to indicate the ability of the 
computer to calculate accurate values of the S, C, and F functions 
throughout the range of axial load defined by the parameter @% = ‘\pr? /er 
For a standard precision of eight significant digits in all calculat- 
ions, very inconsistent results were evident in all functions for val- 
ues of % less than about twelve degrees, Knowing that at 9 = 0, 

S = 4.0, C = 0.5, and F = 0,083, and knowing also the calculated values 
of S, C, and F at g = 12°, linear functions were used as replacements 
to the actual functions in the range 0<f<12°, The use of simple linear 
functions is justified since the actual functions exhibit such small 
percentage variation in this range (see figures 2-2,2-3,2-4) as com- 


pared to the remaining portions of the curves, 


Before deciding on the use of the standard precision of eight 
significant digits for the running of the bulk of the programs, com- 
panion programs were run using a precision of twelve digits. Inaccur- 
acies were found only in the fifth significant digit, this being con- 
sidered trivial and unworthy of a refinement in precision past eight 


significant digits. 


The programs proper are contained on the following pages. In- 
put format is described on statement number 2 of each program. To 


run more than one data card, insert a ..I card and a «branch OO40R 
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card before each successive data card, For the convenience of those 
wishing to investigate the details of these programs, or to modify 
them in any way, a list of Fortran II symbols used is given following 


the programs. 
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eel ELASTIC STABILITY OF FRAMES 
ool FRAME WITH HINGED BASFS — SMALL DEFLECTION THEORY 
eeLOAD FOFTRAN EXECUTE 
READ 2sHsREI sRPsAND 
2 FORMAT ( F6e2sF6025Fb0e25F be?) 
G=H/RETI 
PUNCH 4 
4 FORMAT(38H1 ELASTIC STABILITY OF FRAMES PART 1R) 
PUNCH 2sHsGseRPsRFT 
3 FORMAT(////92Xs7HLANDA= 9F50295Xs 7HGAMMA= 9F5e2s5Xs4HRP= oF 6035 
15Xs5HRET= 9F6,2) 
DIMENSION P(50) 
P(1)=0.0 
E=2.718281828459 
PY=3.141592653590 
J=2 
ADDX=ADDN 
TH=ANDX 
1 SINTH=SINF(TH) 
COSTH=COSF(TH) 
SAB=4,0 
SBC=4,0 
SCD=4,0 
CAB=0.5 
CBC=0.5 
CCD#£065 
AA=2e*¥G+4e 
AB=2.6 
AC=2.6 
AD=2-6 
AE=AA 
AF=f7 
AG=-( (PY*##2)*G*H)/126 
AH=AG*(-1.) 
AlL=(PY*¥*2)*(RP+TH) 
YA=2e*G*TH 
YB=YA 
YC=€e*TH 
D=AA*AE#*#A I+AD*#AH*¥AC+AB *¥AF XAG~AC¥AE*#AG—-AB¥AD¥A [-AF#AH¥AA 
THB= (YA*X#AE*A 1+ AD#AH*¥YC+YB*¥AFXAG—-YC#AFXAG~YR*AD*AI—-AFHAH#YA) /D 
THCH=(AAFYB¥AI+YA#AH#AC4+ABHYC#HAG—-AC#Y BHAG—-ARHYAHAT—-YC#AH#AA) /ND 
Py =(AA*#AE#YCH+AD#FYBR¥AC4+ABFAFHYA-AC#HAE #YA—-ABHADHYC-AFHYR¥AA)/D 
19 RAB= P1/26-(SRCO#(1.4+C8C)¥( THR4+THC) )/( (PY#*2) *¥G*H) 
PAB=PY*SQRTFI(ARSF (RAB) ) 
IF (RAB)10091019101 
100 IF (PAB-(12¢*PY)/1802)14914912 
14 SAB=4,04+(PAB¥(180.)*(0.0039965) )/(12e*PY) 
CAB=0,5-( PAB¥(180.)* (0200110865) )/(12e#PY) 
GO TO 20 
12 EP=E**PAB 
STHPAB=0.5*(EP—-(1-e/EP)) 
COHPAB=0e5*(EP+(1-2/EP)) 
SAR=(PAR*¥STIHPAR—(PAR*¥#2) *COHPAR) /(2¢*COHPAB-2.—PAR¥SIHPAR) 
CAB=(PAB-SIHPAB) /(SIHPAB-PAB*COHPAB ) 
GO To 20 
101 IF(PAB—(12e#PY)/18006)15915912 
15 SAB=4,.0-(PAR¥(180.)*(020060669) )/(12e*PY) 
CAB=0.54+(PAB#(180.)* (02600110366) )/(12e*PY) 
GO TO 20 
13 SINPAB=SINF(PAB) 
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20 


200 
24 


22 


201 
25 


30 


300 
34 


32 


301 
35 


40 


COSPAB=COSF (PAB) 

SAB=(PAB*¥SINPAB-(PAB*¥*2)*COSPAR) o-2eX* pa 

SUTRAG=SINPAB) /(5INPAB=PAB¥COSPAG, =A er eee 

GO TO 20 

RBC=-(P1*G#H/26+SBC#(1 6+C8C) *( THR+THC) / (PY*#2))*TH 
—(SCD*G#H* (1 --CCND**2)*( THC-TH)/(PY*#2)) 

PBC=PY*SORTF(ASSF(RRC)) 

IF (RBC) 20092013201 

IF (PBC-(120e*®PY)/1800)24924 522 

SBC=4e90+(PBC#(180.)*(0.0039965) )/(12,.*PY) 

CBC#065-(PBC¥(180.)*(0.00110865) )/(12e*PY) 

FBC=0.08333333-(PBC* (180. )*(0.00034900))/(126*PY) 

GO TO 30 

FP=E*#*PBC 

EP2=Ex*(PBC/2,.) 

SIHPBC=0.5*(FP-—(1./FP)) 

COHF BC=0-5*(FP+(1./FP)) 

TAHPBC=(EP2—(146/FP2))/(FP24+(1./5P2)) 

SBC=(PBC*#¥SIHPRC-(PBRC##*2 ) ¥COHPRC) /(2.*COHPRC-2.—-PRC¥STHPRC) 

CBC=(PRC-SIHPRC)/(STHPRC-PRC¥COHPRC) 

FBC=(1le—-(2e*TAHPBC)/PBC)/(PBC##2) 

GO TO 30 

IF (PBC—(12¢*PY)/180.)25925923 

SBC=4.0-(PBC¥(180.)*(0.0060669))/(12-*PY) 

CBC=0.,5+(PBC#(180.)*(0.00110366) )/(12e*PY) 

FBC=0.08333333+(PBC#(1806) *(02.00036456))/(12e*PY) 

GO TO 30 

SINPBC=SINF(PPC) 

COSPBC=COSF (PBC) 

TANPBC=SINF(PRC/2.)/COSF(PBC/2,.) 

SBC=(PBRC¥SINPRC-—(PRO¥*#2)*COSPRCO)/(2.-26*COSPRC-PRCXSINPRC ) 

CBC=(PBC-SINPRC)/(SINPRC-PBC¥COSPRC) 

FBC=((2e*TANPBC)/PRC-1.)/(PBC#*2) 

GO TO 30 

RCD= P1/20t+(SBC#(16+CRC)*( THBR+THC) )/((PY**2) ¥G#H) 

PCD=PY*SQRTF(ABSFIRCD)) 

IF (RCD) 30053019301 

IF (PCD—(12e*PY)/180¢) 24934932 

SCD=4004(PCD¥(18006)* (020039965) )/(12e*PY) 

CCOD=0.5—-(PCD#¥(180.)%*(90.00110865))/(12e*PY) 

GO TO 40 

FP=E#*PCD 

STHPCD=0.5*(EP-(1./FP)) 

COHPTD=0.5*(FP4+(12/7EP)) 

SCD=(PCD*SIHPCD—(PCD*¥*2) *COHPCD) /(2¢*COHPCD-2¢-PCD*SIHPCD) 

CCD=(PCD-SIHPCD) /(SIHPCD-PCD*COHPCD) 

GO To 40 

IF (PCD—(12e*PY)/180¢) 35935933 

SCD=400-(PCD¥(1800)* (020060669) )/(12e*PY) 

CCOD=0054(PCD¥(1800) * (0000110366) )/(126*PY) 

GO TO 40 

SINPCD=SINF(PCD) 

COSPCD=COSFI(PCD) 

SCDe (PCD¥SINPCD—(PCD##2) ¥COSPCD) J (26-20 *COSPCD—PCD¥SINPCN) 

CCD=(PCHN-—SINPCD) /(SINPCD-—PCD*COSPCD) 


GO TO 40 
AA=SAB¥G¥*(1e-—(CAB¥¥2) )4+SBC 
AB=SBC#CBC 
AC=SAB¥(1e—-(CAB**2) ) 

AD=AB 
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4 FORMAT(38H1 ELASTIC STABILITY OF FRAMES PART 1A) 
PUNCH 39H sGsRPsREI 
3 FORMAT(////32Xs7HLANDA= 9F5e¢295Xs7THGAMMA= 9F5e295Xs4HRP= oF 603 
15Xs5HREI= »F6e2) 
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GO TO 20 

SINPAB=SINF (PAR) 

COSPAB=COSF (PAR) 

SAB=(PAB¥SINPAR—(PAB¥*2) *COSPAR) / (26-26 *COSPAR-—PAR¥SINPAR) 

CAB=(PAB-SINPAB)/(SINPAS8—PAR¥COSPAR) 

GO TO 20 

RBC=—(P1*¥G*H/2.4+SRC#(1.+CB8C) *( THBR+THC) /(PY**2))¥*TANTH 
—(SCD*G#H*(1.—-CCD*¥*2)*( THC-TANTH) /(PY**2))/COSTH 

PBC=PY*SORTFI(ABSFI(RBC)) 

IF (PBC) 20092013201 

IF (PBC-(12e#PY)/1800)24924522 

SBCH#4e04+(PBC#(180.)*(0.0039965))/(12—_*PY) 

CBC=0*5-—(PBC¥(180.)*(0.00110865))/(120e*PY) 

FBRC=0.08233332-—(PBRC#(180~6)#(0.00034900))/(12e*PY) 

GO TO 30 

EP=E*x*#PBC 

EP2=E#*(PBC/2,) 

STHPBC=0.5*(EP-(1./FP)) 

COHPBC=0-5*(EP4+(1./EP) ) 

TAHPBC=(EP2-(1-/EP2))/(EP24+(1./EP2)) 

SBC=(PBC*¥STIHPBC-(PBC##*2) *COHPBC) /(26*COHPRC—-2.—-PRC¥SIHPRC) 

CBC=(PBC-SIHPBC)/(STHPBC-PRC¥COHPBC ) 

FBC=(1e-(2e*TAHPBC)/PRC)/(PBRC#¥2) 

GO TO 30 

IF (PBC—(12¢*PY)/1800)25925923 

SBC=4-0-(PBC#(180.)*(0.0060669) )/(12e*PY) 

CROH06.54+(PRCO¥(1806.)*(0.00110366))/(126#PY) 

FRC=0,.0833332334+(PRC#(1806.)*(0.00036456))/(12e*PY) 

GO TO 30 

SINPBC=SINF(PBC) 

COSPBC=COSF (PRC) 

TANPBC=SINF(PBRC/2.)/COSFIPBC/2.-) 

SBC=(PBC¥SINPBC-(PBC#¥*2) *COSPRC)/(26-26*COSPBC-PBRC¥SINPRC) 

CBC=(PBC-SINPBC)/(SINPRC-PBC*¥COSPBC) 

FBC=((2e*TANPBC) /PBC—-1e)/(PBC#¥*2) 

GO TO 30 

RCD=(P1/20t+(SRC#(104+C8C) *( THR+THC) )/( (PY*¥*2)*G*H) )/COSTH 
+(SCD¥(1e-CCD**2) *( THC-TANTH) *TANTH) / (PY**2 ) 

PCD=PY*SORTFIARSFI(RCD)) 

IF (RCD) 32009301 9301 

IF (PCD-(12e*PY)/1806) 34934932 

SCD=4004+(PCD¥(180.)* (020039965) )/(12e*PY) 

CCN=0,5-(PCH¥(1806.)*(0.00110865))/(12e*PY) 

GO TO 40 
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SCD=(PCD*¥STHPCD—(PCD*¥*2) *COHPCD) /(2¢*COHPCD—-26-PCD*SIHPCD) 

CCD=(PCD-SIHPCD)/(SIHPCD-PCD*COHPCD) 

GO TO 40 

IF (PCD—(120#PY)/1806) 325935333 

SCD=4-,0-( PCD¥(180.)* (020060669) )/(126*PY) 

COD=I0954+(PCD¥( 1800 )* (0200110366) )/(120*PY) 

GO TO 40 

SINPCD=SINF(PCD) 

COSPCD=COSF(PCD) 

SCD=(PCD¥SINPCD—( PCD¥#*2)*COSPCD)/(2¢-2e¢*COSPCD-PCD¥SINPCD) 

CCD=(PCD-SINPCD)/(SINPCD-PCD*COSPCD) 

GO TO 40 

AA=SAB¥G¥(1e—(CAB¥*2) )+SBC 
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AF=SCD¥(1.-(CCD¥*2)) 
AG=-(SBC#(1.-CBC) *G#¥H*(PY#*2)#FRBC)/2, 
AH=AG*(—-1.) 
AT=(PY*¥*2)*(RP¥COSTH4+SINTH) 
YA=SAB*#G¥(1.—-(CAB¥*2) )*¥TANTH 
YB=SCO*G¥(1.—-(CCD#¥*2))*TANTH 
YC=(SAB*¥(1e—(CAB**2 ) )4+SCD¥(1 e-—(CCD**2)))*TANTH 
D=AA*AE #A IT+AD#AH¥FAC+AR RAF RXAG—-ACHAEXAG—-AB¥AD#AIT—AFXAHFAA 
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THCH(AA¥YBH¥AIT+YA¥AHXACTH4ABRKYCHAG-AC#YR¥XAG~-ARKYA#KAI-YC#HAH¥AA) /D 
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PPP=(ABSF{P2-P1))/P1-0.005 
IF (PPP)51951%50 
50 P1=P2 
GO TO 10 
51 P(J)=P2 
THD=(164+CCD) #*TANTH-CCD*¥THC 
IF (ABSF (P2)-0.6)49949970 
49 PUNCH 52sTHsP2sTHBsTHCsTHD 
52 FORMAT(// 2X9 THTHETA= 9Flle8s9H ALPHA= sF1l00e627H THR= sF6039 
17H THC= 9F60e397H THD= 9F 603) 
PUNCH 53sRABsRRC eRCD 
53 FORMAT(41Xs5HRAR= sF60e397H RBC= oF 6e397H ROCD= 9F 603) 
PPPP=P(J)-P(J-1) 
IF (PPPP)70970960 
60 IF (J-5)69963964 
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GO TO 69 
64 IF (4J-13)69265969 
65 ADDX=3«4*ADD 
GO TO 69 
69 J=J+1 
TH=TH+ADDX 
IF (TH-0 298) 88588870 
88 GO TO 1 
70 CALL EXIT 
62 END 
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eel ELASTIC STABILITY OF FRAMES 
eel FRAME WITH FIXED BASES — SMALL DEFLECTION THEORY 
eeLOAD FORTRAN EXECUTE 
READ 2sHsREI sRPs AND 
2 FORMAT ( F6e29F6e29Fb02 F602) 
G=H/RET 
PUNCH 4 
4 FORMAT(38H1 ELASTIC STABILITY OF FRAMES PART 2B) 
PUNCH 39HsGsRPsPRF] 
3 FORMAT(////92Xs7HLANDA= 9F50235X%s 7HGAMMA= 9F5e2e5Xs4HRP= gF6039 
15XsSHRET= 9F6.2) 
DIMENSTON P(50) 
P(1)=0.0 
E=2.718281828459 
PY=322141592653590 
J=2 
ADDX=ADD 
TH=ADDX 
1 SINTH=SINF(TH) 
COSTH=COSF(TH) 
SAB=4,0 
SBC=4,0 
SCD=4.0 
CAB=0,5 
cBc=0.5 
CCD=965 
AA=4e%*G+4. 
AB=2. 
AC=66 
AD=2.6 
AE=AA 
AF=AC 
AG=—( (PY##2)#G*H)/12, 
AH=4£G*(-1.) 
AT=(PY*#*2)*(RP4+TH) 
YA=€ «*G*TH 
YB=YA 
YC=24e*TH 
D=AA*AE*A IF+AD#AH¥AC+ABXAF #AG—-AC#AE XA G—AB*AD¥A I—-AF#AH*#AA 
THB=(YA*AE*AI+AD*¥AH* YC4+YB*¥AF*XAG—-YC#AFRAG—-YBXADHAI-AFHXAH¥YA) /D 
THCH(AA*¥YB#AI+YA*#AH*¥AC+AB*YC#AG—-AC#YBR¥AG-AB¥YA*#A IT-~-YC¥AH¥AA)/D 
PT =(AA*AE*XYC+AD*YB*ACH+AR *AFRYA—-AC#AE ¥YA—-AB¥AD#HYC-AF#XYR¥AA)/D 
10 RAB= P1/2e-(SBC#¥(1e+CRC)*( THB+THC) )/(( PY¥*2)*G#*H) 
PAB=PY*SOQRTF(ABSFI(RAB)) 
IF (RAB)10091019101 
100 TF(PAB—(120¢*PY)/1804.)14914912 
14 SAB=4.0+4+(PAB*¥(180.)*(0.0039965) )/(12e*PY) 
CABH0.5—-(PAB¥(180.)*(0200110865) )/(12e*PY) 
GO TO 20 
12 EP=E**PAB 
SITHPAB=0.5*(EP-(1-/EP)) 
COHPAB=0.5*(EP+(1./EP)) 
SAB=(PAB*SIHPAB-(PAB##2) *COHPAB) /(2¢*COHPAB—-2e—-PAB¥STHPAB) 
CAB=(PARB-SIHPAR) /(SITHPAR-PAB*COHPAB ) 
GO TO 20 
101 TF(PAB—(12e#PY)/1804)15915913 
15 SAB=4,0-(PAR¥(180,.)*(0.0060669))/(12—*#PY) 
CABH=0.54+{(PAR¥(180,)*(0,00110366) )/(12e*PY) 
GO TO 20 
13 SINPAB=SINF(PAB) 
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COSPAB=COSF(PAB) 

SAB=(PAB*SINPAB-—(PAB¥*2) *COSPAB) /(2.-26*COSPAB—-PAB*¥SINPAB) 

CAB=(PAB-SINPAB)/(SINPAB-PAB*¥COSPAB) 

GO TO 20 

RBC=—(P1*G*H/26+SBC#¥(1.+C8C)*( THR+THC) /(PY*¥*2))*TH 
—( SCD*¥G¥H¥(16+CCD) *( THC-2.*TH)/(PY**2) ) 

PBC=PY*SQRTF(ARSF(RBC) ) 

IF (RBC) 20092019201 

IF (FP 8C-(120*PY)/180.)24924922 

SBC=4504+(PRC#(180.)*(0.0039965))/(12—-*PY) 

CBC=945—(PBC#(180~.)*(0.00110865) )/(126*PY) 

FBC=0.08333333-(PBC#*(1800)*(0.00034900) )/(12e*PY) 

GO TO 30 

EP=E**#PBC 

EP2=E*#*(PBC/2.) 

STHPBC=0.5*(EP-(1./EP)) 

COHPBC=0.5*(EP+(1-e/EP)) 

TAHPBC=(EP2-(16/FP2))/(EP24+(1./FP2)) 

SBC=(PBC¥SIHPBC—(PBC*¥*#2) *COHPRC) /(2.*COHPRC-26—-PRC¥STHPTC) 

CRC=(PBC-SIHPRC)/(STHPRC-PRC*¥COHPRC ) 

FRC=(1e—-(2e*#TAHPSC)/PRCO)/S(PRC##2) 

GO To 30 

IF (PBC—(120¢*PY)/180¢0)25925323 

SBC=4,0-(PBC*¥(180.)*(0,0060669))/(126*PY) 

CBCH0554+(PBC#(1806.)*(0.00110366) )/(12e*PY) 

FBC=0.0833232334+(PBC¥(1800)*(94.000326456))/(12e*PY) 

GO TO 30 

SINPBC=SINF(PBC) 

COSPBC=COSF(PBC) 

TANPBC=SINF(PBC/2.)/COSFI(PBC/2.) 

SBC=(PRC#¥SINPBRC—(PRO¥*2)*COSPRC) /(2.-26*COSPBRC—PRC¥SINPRC ) 

CRC=(PRC-SINPRC)/(SINPRC-PRC¥COSPRC) 

FBC=((2e¢*TANPBC)/PRC-12.)/(PRC#*2) 

GO TO 30 

RCD= P1/2e+(SBC#¥(16+C8C) ¥*( THB+THO) )/( (PY**2 )*#G*H) 

PCD=PY*SQRTF(ABSFI(RCD)) 

IF (RCD) 30093015301 

IF(PCD—(12e#PY)/1806¢) 34934532 

SCD=4604+(PCD¥(1800)* (020039965) )/(12e*PY) 

CCD=065-—(PCD¥(1806¢)*(0.00110865) )/(12e*#PY) 

GO TO 40 

EP=F##PCD 

STHPCD=0.5*(EP-(1./EP)) 

COHPCD=0.5*(FP+(1./FP)) 

SCD=(PCD¥SIHPCD—(PCD¥¥2) *COHPCD) / (2. *COHPCD-2.—-PCD¥SIHPCD) 

CCD=(PCD-—SIHPCD) /(SITHPCD-—PCD*COHPCD) 

GO To 40 

IF (PCD-(120#PY)/180¢) 35935933 

SCD=440-( PCD¥ (1806) *(020060669))/(12e*PY) 

CCOD=£06 5+(PCD¥(1806)*(02.00110366) )/(126*PY) 

GO TO 40 

SINPCD=SINF(PCD) 

COSPCD=COSFI(PCD) 

SCD=(PCD*SINPCD—(PCD#*¥2)*COSPCD) /(2e-2e*COSPCD—PCD¥SINPCD) 

CCD=(PCD—SINPCD) /(STNPCD-PCD¥COSPCD) 

GO TO 40 

AA=SAB¥G+SBC 

AB=SBC¥CBC 

AC=SAB¥(1¢+CAB) 

AD=AB 
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AE=SCD*G+SBC 

AF=SCD*(1.2+CCD) 

AG=-(SBC¥*¥(1.-C8C) *G*¥H*¥(PY*#*2)*FBC)/2. 

AH=AG¥(-1.) 

AT=(PY*#*2)*(RP+TH) 

YA=SAB#¥G*(1-6+CAB)*TH 

YR=SCD¥G¥(1.+C00CND)*TH 

YC=(2.*SAB¥(1.¢4+CAB)+2.*SCN¥(1.2400N) )*TH 

D=AA*AE *#A ItFAD*#AH¥AC+AB #AF XAG—-ACHAE*XAG-AB#¥AD¥AI-—AFXAH*FAA 
THBR=(YA#AE*# A T+AD*#AH*¥YC4+YS*AF ¥AG-YC#HAFXAG—-YR¥AD#A IT-AFX¥AH*¥YA)/D 
THC=(AA*YB¥AI+YA*#AH*FAC4+ABRKYC#AG-AC#YR¥AG-AB¥YA*AI-YC#AH#AA) /D 
P2=(AA* AEX YC+AD*YB*¥AC4+AB RAF HY A-—AC#KAE#YA-AB¥AD¥#YC-AFHYB*AA) /D 
PPP=(ABSF(P2-P1))/P1-0,005 

IF (PPP)51951950 

P1=P2 

GO TO 10 

P(J)=P2 

PUNCH 52sTHsP2sTHBsTHC 

FORMAT(//5s2X9 THTHETA= 9F1l1le8s9H ALPHA= sF100e6s7H THB= »sF6e3>9 


17H THC= »F603) 


PUNCH 53sRABsRRCeRCD 
FORMAT(41Xs5HRAB= oF60397H RBC= sF603297H RCD= oF 603) 
PPPF=P(J)-P(J-1) 

IF (PPPP)70970%60 

IF (J-5)69 963964 
ADDX=246*ADD 

GO TO 69 

IF (J-13)69265969 
ADDX=34*ADD 

GO TO 69 

J=J+1 

TH=TH+ADDX 

IF (TH-1205) 885988970 
GO TO 1 

CALL EXIT 

END 
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eel ELASTIC STABILITY OF FRAMES 
eel FRAME WITH FIXED BASES — LARGE DFFLECTION THEORY 
eeLOAD FORTRAN EXECUTE 
READ 2sHsREI »RP»s AND 
2 FORMAT (F6 e2sF6e29FGe29Fbe?) 
G=H/RET 
PUNCH 4 
4 FORMAT(38H1 ELASTIC STABILITY OF FRAMES PART 2A) 
PUNCH 39HsGeRPosREI 
3 FORMAT(////52Xs7THLANDA= 9F50295X97HGAMMA= 9F50295X94HRP= 
15Xs5HREI= F602) 
DIMENSTON P(50) 
P(1)=0e0 
EF=22718281828459 
PY=3.141592652590 
J=2 
ADDX=ADD 
TH=ADNX 
1 SINTH=SINF(TH) 
COSTH=COSF(TH) 
TANTH=SINTH/COSTH 
SAB=4.0 
SBC=4.0 
SCD=4.0 
CAB=0.5 
CBC=0.5 
CCD=045 
AA=4e*G+4. 
AB=2e6 
AC=6-6 
AD=26 
AE=AA 
AF=fAc 
AG=—((PY*##2)*¥G*H)/12. 
AH=AG¥*(-1.) 
AIT=(PY*¥*2)*(RP*¥COSTH+SINTH) 
YA=6e*G*TANTH 
YB=YA 
YC=24.*TANTH 
D=AA*AE*A I+AD#AHFAC+AB HAF XAG—ACHAEXAG—AB¥AD¥AI—AFXAH*AA 


sF6e39 


THB=l(YA*AE*#AIT+AD#AH#YC+YBRHAF¥XAG—-YCHAFRAG-YB RADHA T—AFHAH¥YA) /D 
THCH=(AA¥*YB¥AT+YA*XAH*¥AC4+AB RY C#AG—AC#Y BEAG=ABRHYA#A T-~YC#HAH¥AA ) /D 
P1=(AA*AEX#YC+AD¥*¥YR¥AC+ABRAFHYA-ACHAE #YA~ABHEADHYC-AFHYRHAA)/D 


10 RAB=(P1/26—-(SRCO#¥(1e+CRC)*( THR+THC) )/ 0 ( PY**2)*G¥H))/COSTH 


1 +(SAB¥(164+CA8)*( THB-2e*TANTH) *#TANTH/ (PY##2) ) 
PAB=PY*SQRTF(ABSF(RAB)) 
IF (RAB)10091919101 
100 IF (PAB-(12¢*#PY)/1800)14914912 
14 SAB=4.04+(PAB¥(180.)*(0.0039965))/(12e*PY) 
CAB=025-(PAB¥(1806)*(0.00110865))/(12e*PY) 
GO TO 20 
12 EP=E*#PAB 
STHPAB=0.5*(EP-(1-2/EP)) 
COHPAB=0,.5*(FP+(1./EP)) 


SAB=(PAB*STHPAR-(PAB**2) *COHPAB) /(2.*COHPAB-2.—-PAR*SIHPAR) 


CAR=(PAR-SIHPAR) /(STHPARB-PAR*COHPAB ) 
GO TO 20 

101 IF(PAB—(12e*PY)/1806)15915913 

15 SAB=4,0-(PAB¥(180.)*(0.0060669))/(126*PY) 
CAB=0.5+(PAB¥(180.)*(02.00110366) )/(12e*PY) 


Del 5 
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GO TO 20 
SINPAB=SINF(PAB) 
COSPAB=COSF (PAR) 
SAB=(PAB*SINPAB—(PABR¥*¥2) ¥COSPAR) /(2.-2e*COSPAR—PAR¥SINPAR) 
CAB=(PAB-SINPARB)/(SINPAR—-PAB¥COSPAR) 
GO TD 20 
RBC=—-(P1*G*H/264+SRBC¥(1.4+C8C)*( THR+THC) /(PY*#*2) )*¥TANTH 
— (SCD*¥G*H*¥(1.64+CCD) *( THC-2e*TANTH) /(PY**2) )/COSTH 
PBC=PY*SQRTF(ABSFI(RBC)) 
IF (RBC) 20092019201 
IF (PBC-(12e*PY)/180¢4)24924922 
SBC=4-.04+(PBC¥(180.)*(0.0039965) )/(126*PY) 
CBC=05.5-(PBC¥(180.)*(0.00110865))/(12e*PY) 
FBC=0408333333-(PBC#(1R806¢ )*¥(0200034900))/(12e*PY) 
GO TO 30 
EP=E*##*#PRBC 
EP2=F**(PBC/2,.) 
STHPBC=0.5*(FP-(1./FP)) 
COHPBC=0.5*(EP+(1-./FP)) 
TAHPBC=(EP2—(12/FP2))/(EP24+(1./FP2)) 
SBC=(PBC¥SIHPBC-(PRC**2) *COHPRC) /(26*COHPBC-2-.-PRC¥SITHPRC) 
CBC=(PBC-SIHPBC)/(SIHPBC-PBC*COHPBC) 
FBC=(1e-(2e*TAHPBC) /PBC)/(PBC##2) 
GO TO 30 
IF (PBC—(12e¢#PY)/1802)25925923 
SBC=4.0-(PBC*(180.)*(0.0060669) )/(126*PY) 
CBCH=0,54+(PBRC¥(180,)*(0.00110366) )/(12e*PY) 
FBC=0.0833333234(PRC#(1R8O6)*(0.00036456) )/(12e*PY) 
GO TO 30 
SINPBC=SINF(P8C) 
COSPBC=COSF (PBC) 
TANPRC=SINF(PRC/2.)/COSFI(PRC/2.) 
SRCH(PBC¥SINPRC—-(PBRC#¥#2)*COSPBC)/(26-26*COSPBC—-PRC¥SINPRC) 
CBC=(PBC-SINPRBC)/(SINPRC-PBC*COSPBC) 
FRC=((2e* TANPBC) /PBC-1.)/(PBC##2) 
GO TO 30 
RCD=(P1/204+(SBC¥(164+CRC) *( THR+THC) )/((PY**2)*G#H) )/COSTH 
+(SCD¥(1e¢+CCD) *( THC—-2e*TANTH) #TANTH/ (PY**2) ) 
PCD=PY*SQRTFI(ABSFI(RCD)) 
IF (RCD) 320093015301 
IF (FPCD—(120e*PY)/180 6) 324934932 
SCD=4e)04( PCDN¥(1806.)* (060039965) )/(12e*PY) 
CON=0.5—-(PCD¥(180-.)* (0200110865) )/(12e*PY) 
GO TO 40 
EP=F— **PCD 
SIHPCD=0.5*(EP—-(1./EP)) 
COHFE CD=0-e5*(EP+(1-¢/EP)) 
SCD=(PCD¥SIHPCD-(PCD**2)*COHPCD) /(2¢*COHPCD-2-e-PCD*SIHPCD) 
CCD=(PCD—SIHPCD) /(SIHPCD-PCD*COHPCD) 
GO TO 40 
IF (PCD-—(12e#PY)/1806) 25935933 
SCD=45,0-( PCD¥(1806.)* (0209060669) )/(12e*PY) 
CONH=0654( PCDO¥(180~6)* (0200110366) )/(12e*PY) 
GO TO 49 
SINPCD=SINF(PCD) 
COSPCN=COSFI(PCD) 
SCD=(PCD¥SINPCD—(PCD¥*¥2)*COSPCD)/(2e-2e*COSPCD—PCD*SINPCD) 
CCD=(PCD-SINPCD)/(SINPCD-—PCD*COSPCD) 
GO TO 40 
AA=SAB#*G+SBC 
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AB=SBC*¥CBC 

AC=SAB¥(1¢+CAB) 

AD=AB 

AE=SCD*¥G+SB8C 

AF=SCND#¥(12+CCD) 

AG=—-(SRC#¥(1.-CRC) #G#H*(PYX*#2)#FRO)/2, 

AH=AG*(-1.) 

AT=(2Y##2)*(RP¥COSTH+SINTH) 

YA=SAB*G*(1¢+CAB)*TANTH 

YB=SCD*G#(1.2+CCD)*TANTH 
YC=(2e*SAB*¥(164+CAB)+2e*SCD¥(1.+CCD) ) *TANTH 

D=AA*AE*A I +AD#AH*¥AC+AB AF XAG—-AC#AE*XAG—AB#RAD*A I —AFXAHFAA 
THB=(YA*#AE#A I+ AD*AH*YC4+YBXAF#XAG-YC#AFE#AG—-YBXAD¥A I-AF¥AH#YA) /D 
THC=(AA*#YBHAT+YA#AH*AC4+AB#¥YC#AG—-AC#YPR#¥AG—-ABHYA#AI—YC#AH#AA) /D 
P2=(AAXAEX# YC+AD#YB*#AC+ABHAF HY A—-AC#AE #YA—-AB¥AD¥YC-AF#HYBEAA) /D 
PPP=(ABSF(P2-P1))/P1-0.005 

IF (PPP)51251950 

Pl=P2 

GO TO 10 

P{J)=P2 

IF (ABSF (P2)-244)49949370 

PUNCH 52sTHsP2sTHB»s THC 

FORMAT(//s2X9 THTHETA= 9F11e8s9H ALPHA= 9F10e6s7H THB= »F6e39 


17H THC= F603) 


PUNCH 52sRABsRBCe»RCD 
FORMAT(41XsS5HRAB= sF6e397H RBC= sF603297H RCD= F603) 
PPPP=P(J)-P(J=-1) 

IF (PPPP)70970960 

IF (J-5)69 963964 
ADDX=24*ADD 

GO TO 69 

IF (J-13)69965969 
ADDX=3e*ADD 

GO TO 69 

J=J+1 

TH=TH+ANDX 

IF ( TH-1205)88288970 
GO TO 1 

CALL EXIT 

END 
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ELASTIC STABILITY OF FRAMES 
FRAK.E WITH PARTIAL BASE FIXITY 
eeLOAD FORTRAN EXECUTE 


READ 2sHeREI eR Ps ADDF 
FORMAT (F6629FSe2sF6029F6 
G=H/REI 

PUNCH 4 


FORMAT(328H1 ELASTIC STABILITY OF FRAMES 


PUNCH 3sHsGsRPsRFI5F 
FORMAT (////32Xs7HLANDA= 


DIMENSION P(50) 
P{(1)=0,0 
E=2.718281828459 
PY=3.141592653590 
J=2 

ADDX=ADD 

TH=ADDX 
SINTH=SINF(TH) 
COSTH=COSF(TH) 
SAB=4.0 

SBC=4,.0 

SCD=4,0 

CAB=0,.5 

CBC=0,.5 

CCD=0.5 
AA=—-60*G-8 eo *F—-Re-1lLOe *F/G 
AB=—4e-4e*F-8e%*F/G 
AC=-5e-l18 e*F/G 
AD=AB 

AE=AA 

AF=AC 
AG=—((PY#*2)*G*H)/12.6 
AH=AG¥(-1.) 
AT=(PY#*2)*(RP+TH) 
YA=(-6e*G-18e)*TH 
YB=YA 

YC=-12e*TH 


9F60293Xs5HFETA= 


o29Fbe2) 


9F5e293X 
9F6e2) 


» THGAMMA= 


(SMALL DEFLECTION THEORY) 


PART 3R) 


9F5e293X94HRPH 


D=Afh FAE XA T+AD¥*¥AH*AC+AB RAF XAG—-AC#AE XA G—-AB*¥AD*AI-—AFXAH*AA 
THA=(YA*#AE*#A I1+AD*¥AH*¥YC+YBXAF*AG—-YC#AF*XAG—-YB*AD*A I-AF*XAH*YA) /D 
THD= (AA*#YB*¥AIT+YAXAH*AC+ABRYC#AG—-AC#Y B¥AG—-AB#YA*A I-YC#AH#AA)/D 
PLT=(AA*®AE*#YC+AD*#YBXAC4+ABKAF HY A-ACHAE HY A-AB*XAD¥YC-AFHYBRAA)/D 
THR=3 6 *TH—( Ge #GH+4 0 #F )*THAS (20 *G)—2e*FRXTHD/ (20*G) 
THCH=B oF TH—( 4a *¥GH+4 Oe FF) FTHD/ (20 GG) —2eFFRTHAS(20%*G) 
RAB= P1/2e-(SBO#¥(164+CRC)*( THR4+THC) )/( (PY*®*2 )¥G*H) 


PAB=PY#*SQRTF(ARSF(RAB)) 
IF (RAB)10091019101 


IF (PAB-—(12e*PY)/1802.)14914912 
SAB=4.0+(PAB*¥(180.)*(0.0039965) )/(12e*PY) 
CAB=0.5-(PAB*¥(180.)* (0200110865) )/(12e*PY) 


GO TO 20 
EP=E*#*PAB 
STHPAB=0e5*(EP-(1.2/EP)) 
COHPAB=0.5*(EP+(1./EP)) 


9F5e3°9 


SAB=(PAB*SIHPAB-— (PAB **2 ) *COHPAB) /(2¢*COHPAB—2.-PAB*¥SIHPAB) 
CAB=(PAB-SIHPAB) /(SIHPAB-PAB*COHPAB ) 


GO TO 20 


IF (PAB—(12¢#PY)/1806)1591591% 
SAB=4.0-(PAB¥(180.)*(0.0060669) )/(126*PY) 
CAB=0.5+(PAB¥(180.)*(0.00110366) )/(12e*PY) 
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GO TO 20 
SINPAB=SINF(PAP) 


COSPAB=COSFI(PAB) 

SAB=(PAB*¥SINPAB- (PAB¥*2)*COSPAB)/(2,.-2.*COSPAR-PAR¥SINPAB) 

CAB=(PAB-SINPAB)/(SINPAB-PAR¥COSPAR) 

GO TO 20 

RBC=—(P1*G*H/2 6+SBC*(1.6+CRC)*( THR+THC) S/(PY*#*2) )*TH 
— (SCD*G*H*(1.4+CCD) * ( THC+THD-2 0*TH) )/( PY**2 ) 

PBRC=PY*SQRTF(ARSF(RRC)) 

IF (RBC)20092019201 

IF (PRC—(12e*#PY)/180.)24924922 

S§BC=4.04+(PRCO#(180.)*(0.0039965))/(12.*PY) 

CBC=0.5=—(PRC#(180.)*(0.00110865))/(126*PY) 

FBC=0.08333333-(PRC# (180. )*(0.9000324900) )/(126*PY) 

GO TO 30 

FP=E**PBC 

EP2=E¥*(PBC/2.) 

STHPBC=0.5*(EP-(1./EP)) 

COHPBC=0-.5*(EP4+(1./FP)) 

TAHPBC=(EP2—(16e/EP2))/(FP24+(16/EP2)) 

SBC=(PBC¥SIHPBC-(PRC¥*2) *COHPRC) /(2.*COHPRC-26—-PRC¥STHPRC) 

CBC=(PRBC-SIHPRC)/(STHPRC-PBC*COHPRC ) 
FRC=(1e-(2e*TAHPRC)/PRC)/(PRO#¥2) 

GO TO 30 

IF (PBC-(120e#PY)/180.)25925323 

SBC=4.0-—( PRBC¥(180-.)*(0.0060669) )/(126*PY) 

CBC=H=0.54+(PBC*(1806.)*(0.001103646) )/(12e*PY) 
FBC=0,.083333334+(PBC#¥(1800 )*(0.900326456) )/(12e*PY) 

GO TO 30 

SINPBC=SINF(PBC) 

COSPBC=COSF (PBC) 

TANP3C=SINF(PBC/2.)/COSFI(PBC/2.) 

SBC=(PBC*¥SINPBC-(PBC¥*2) ¥COSPBC)/(2.-2e*COSPRC-PRC¥SINPRC) 
CBC=(PRC-SINPRC)/(SINPRBC-PBC¥COSPRC) 
FRC=((2e*TANP8C)/PRC-1-.)/(PBC#¥2) 

GO TO 30 

RCD= P1/26+(SBC#(164+CR COC) *( THR+THO) )/S( (PY*##2 )¥G#H) 
PCD=PY*¥SORTFI(ARSFI(RCD)) 

IF (RCD) 320093013301 

IF (PCD—(120#PY)/1809¢) 24934932 
SCD=4.04+(PCD¥(180.)* (020039965) )/(12e*PY) 
CCDH0.5-—( PCD¥(1806)*(0.00110865) )/(126*PY) 

GO TO 40 

EP=— **¥PCD 

SIHPCD=0.5*(FEP-(1./FP)) 

COHPCD=0.5*(FP+(1./FP) ) 
SCD=(PCD¥SIHPCD—(PCD¥#¥2) *COHPCD) / (26 *COHPCD—-24—-PCD¥STHPCD) 
CCD=(PCN-SIHPCD) /(STHPCD—-PCD*COHPCD) 

GO TO 40 

ITF (PCD-(120*PY)/180¢) 259359232 
SCD=4.0-(PCD¥(180.)*(0.0060669))/(124*PY) 
CCD=0.54+(PCN¥*(180.)* (9600110366) )/(12e*PY) 

GO TO 40 

SINPCD=SINF(PCD) 

COSPCD=COSF(PCD) 
SCD=(PCD¥SINPCD-(PCD##*2)*COSPCD)/ (20-20 *COSPCD—PCD*SINPCD) 
CCD=(PCD-SINPCD) /(SINPCD-PCD*COSPCD) 

GO TO 40 
AA=SAR¥G¥CAB—( SAR¥G4SRC) (SAR #G4+40*F )/(SAR¥G¥CAB) 
—2 oe ¥F*XSBCHCBC/(SCD*¥G*CCD) 
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AB=-SBC¥CBRC¥(SAB¥G4+4.#F) /(SAR¥G¥CABP) 
1 -—(SBC+SCD*G)#*(2e*F)/(SCD¥G¥CCD) 
AC=SAB¥ (1 6+CAB)—(1e+CA8)¥*(SAB¥G+40*F )/(G¥CAB) 
1 —€C1e+CCD)*(2e*F)/(G¥CCD) 
AD==(SAB*G+SBC)¥*(2.*F)/(SAB¥G¥CAB) 
1 -SBC#CBRC¥(SCD¥G4+4e*F)/(SCD#¥G¥CCD) 
AE=SCD*G*¥CCD—(SBC+SCD¥G) ¥*(SCD#¥G+4e*F )/( SCD¥G¥CCD) 
1 -SBC¥CBC#(2e*F)/(SAB¥G¥CAR) 
AF =SCD¥(10t+CCD)—(10+CCD) (SCD¥G4+4G0*F)/(G¥CCD) 
1 —(1e+CAB)*(2,.%F)/(G¥CAR) 
AG=-SBC#¥(1e¢-CBC) ¥G*H*(PY**#2)*#FRO/2. 
AH=AG*¥(=-1.) 
AT=(PY*¥#2)*(RP+TH) 
YA=(SAB¥G*(16+CAB) -—(SAR*¥G+SRC)¥(16+CAB)/CAR 
1 ~SFC¥CBC#¥(1.64+CCD)/CCD)*TH 
YB=(SCD*G¥(10+CCD)—-(SBRC+SCD¥G)¥(164+CCD)/CCD 
1 -SEC¥CBC¥(1¢+CAB)/CAB)*TH 
¥C=(2e*SAB*(164+CAB)4+2 o*SCD¥(164+C0CD)-SAB¥(16¢+CAB)¥(10+CABR)/CAB 
1 -SCD¥(1¢+CCD) *(164+CC09)/CCD)*TH 
D=AA*AE*A I+AD#AH¥ACH+AB¥AF XAG-AC#AE*XAG—-AB¥AD¥A I-AF*XAH¥AA 
THA=l(YA#AE*A I+ AD#AH*FYC+Y 8B *AF¥AG-YCHAFXAG—-YBXAD¥A I-AFXAH¥YA) /D 
THD=(AA*YB¥AI+YA*AH¥AC4+AB#YC#¥AG-AC#YR¥AG—-AR¥YA¥A T-YC#AH#¥AA) /D 
P2=(AA#AEXYC+AD¥*#YR¥AC+AB¥EAFHYA—-ACHAEXYA—-AB¥AD¥YC-AFHYREAA)/D 
THB=(14+CAB) *TH/CAR— (SAB ¥G4+4 0 *F )¥THAS(SAR¥G¥CAB ) 
1 —-2e*F*THD/ (SAR*¥G¥CAB) 
THCH=(104+CCD) ¥TH/COD—(SCDO¥G44 0 *F) XTHD/(SCD¥GECCD ) 
1 -2e*F*THA/(SCN*#G*CCD) 
PPP=(ABSF(P2-P1))/P1-0.005 
IF (PPP)5151950 
50 Pl=P2 
GO TO 10 
51 P(J)=P2 
PUNCH 529THsP2sTHBRsTHCsTHD 
52 FORMAT(//92Xs THTHETA= 9Flle8s9H ALPHA= sF1l00e697H THB= »F603s 
17H THC= sF6e397H THD= »F602) 
PUNCH 532sRAR2R8BC esPCD 
532 FORMAT(4GIXsSHRAB= »F60297H RRC= sF6e%97H ROD= sF 60) 
PPPP=P(J)-—P(J-1) 
IF (PPPP)70970960 
60 IF (JU-5)69963964 
63 ADDX=24*ADD 
GO TO 69 
64 IF (J-13)69965969 
65 ADD¥Y=3-4*ADD 
GO TO 69 
69 J=J+1 
TH=TH+ADDX 
IF (TH=1205)88528970 
88 GO TO 1 
70 CALL EXIT 
62 END 
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Fortran II Symbols 


Mee T = PY 
A =H 9 = TH 
8 =G 0, = THA 
< = P 0, = THB 
RP = RP ®. = THC 
REI = REI 0, = THD 
e=E 
Increment to © values — ADDX 
Initial ADDX value - ADD 
Sap = SAB Dee = SBC S. = SCD Faao= fee 
C,, = CAB Cy. = CBC C,, = CCD 
Pas RAB Pec = RBC %» = RCD 
D,. = PAB Dae = PBC G., = PCD 
Sin 9. SINPAB Cos 9, COSPAB 


= = 
Sin J,, = SINPBC Cos Jac = COSPBC Tan Yee = TANPBC 
= = COSPCD 2 


Sinh J,, = SIHPAB Cosh Y,, = COHPAB 

Sinh Jy, = SIHPBC Cosh Jy, = COHPBC ‘Tanh Jy. = TAHPBC 
Sinh Gp = SIHPCD Cosh G,, = COHPCD 2 

Sin © = SINTH Cos © = COSTH Tan 0 = TANTH 

ef = EP 

e% = EP2 
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